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ABSTRACT 


Material requirements planning is extensively used 
in the industry as the core part of the total manufacturing 
planning and control. To obtain maximum advantage from it, it is 
important that its working be efficient. Many Lot Sizing and 
Scheduling rules, their combinations and their modified forms have 
been tried to give better results in a MRP environment. 

The present work is a simulation study of a MRP 
system to study the effect of coordination between various 
components reaching a assembly centre. In it, it has been 
attempted to reduce the time lag between components so that time 
of waiting at a centre by the components can be saved and by doing 
it the performance of the MRP system can be improved. In order to 
it, a coordination index along with a standard scheduling rule is 
used to decide the priority of components at levels earlier than 
the level at which it gets assembled. 

On application of this approach it was observed 
that it gave significantly better results on important parameters 
like total units lateness and total aggregate lateness. It even 
improved the average capacity utilization of machines. For 
comparison, two other approaches of which one tries to selectively 
reduce lot sizes and the other tries to expedite the late reaching 
components, were considered. 

It is hoped that the developed approach would be of 
some practical help to a practitioner of MRP and would provide a 


new direction for future research work in this area. 
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INTRODUCTION 


In production context, inventory is an idle resource. It 

includes the production material like tools, purchased parts, raw 

\ 

materials, product in process. The fact that the resource is idle 
does not mean it is serving no purpose. It serves as an insurance 
policy against unexpected breakdowns, delays and other 
disturbances that could disrupt the production process. Since cost 
is associated with the disruption of the production process and 
also with maintaining a stock of idle resources, the task of 
operations management is to secure a economic balance between the 
cost of production time loss and the cost of preventing it. 

The most important function of inventory is insulation. 
A reserve of material can be tapped whenever delay in a preceding 
stage threatens to curtail operations in the following stage. 
Material buffers are used to cushion the production process from 
uncertainties of material deliveries, to decouple progressive 
stages of product development from disruption in previous stages, 
and to provide a steady supply of finished output for the unsteady 
demand of customers. 

Operations management, since it involves operational, 
material and financial aspects of a enterprise, is the heart of 
the manufacturing organization. An efficient operations 
management is essential for the health of an enterprise. 

Large manufacturing organizations, such as one engaged 
in aircraft manufacturing, have series of machining and assembly 
operations. Scheduling and inventory control become important 
aspects of operations management in these manufacturing 
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aspects of operations management in these manufacturing 
organizations. 

In an MRP system inventory control and scheduling 
techniques are integrated. It consists of a series of steps which 
start by determining what finished products are needed to meet the 
demand of different time periods, and arrives at a time schedule 
of the finished product components needed at each assembly for 
each time period. 

The material requirements planning system was developed 
to improve the performance of a production system in a complex 
manufacturing setup involving both assembly and manufacturing 
operations. Distinctive characteristics of these complexities make 
the traditional methods of inventory control inadequate. These are, 
(i) the demands and resulting schedules for all components are 
directly dependent on the schedule of requirements for higher 
level assemblies, (ii) Due to limited machine capacities, queues 
are formed before them. (iii) Sequencing decisions affect 
inventory related performance of the shops. The occurrence of 
dependent demands is used by MRP system for coordinating the 
inventory ordering policies for all the components of the 
manufacturing process. Coordination is essential In the case of 
dependent demands as shortage of a single Item can stop the 
assembly process, which leads to low capacity utilization and 
Increase in inventory carrying costs. In MRP environment, the 
demand for the end products as well as intermediate 
products/assemblies is dynamic. 
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1.1 MRP Operation 

A bill of materials describes how a product is made from 
its component parts and assemblies. An MRP system links the 
planned production schedule with the bill of materials needed to 
make the product and examines the manufacturing inventory to see 
which parts and raw materials have to be ordered . By considering 
when various components of the end product are scheduled to be 
produced and the necessary lead times for supply, it time phases 
replenishment orders, so that parts and materials are available 
when they are needed at work stations. Ideally, replenishments 
would arrive at the stations exactly when they are required for 
use, thus minimizing the amount of work in process inventory. 
Since ideals are rarely attainable, the system reexamines the 
production schedule, continuously or periodically, to recognize 
schedule disruptions and irregularities as soon as possible; this 
reduces stockpiles of inventory in the production area. Success of 
adjustments relies on the accuracy of the forecast, bill of 
materials, processing times, routines and inventory records. 

In the total manufacturing control system the position 
of MRP is right in the middle of the total control system. 
Manufacturing control results from forecasts that guide capacity 
planning and inventory control. Master scheduling considers 
available capacity , stock status, forecast demand and customer 
orders to develop a production schedule. The MRP system issues 
planned replenishment orders based on information about job status 
and work in process inventory levels. Production results from 
releasing the planned orders in terms of work loads. Data needed 
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to control operations are obtained by measuring work centre 
output, work in process inventory, and job status. 

The MRP system envelops the entire manufacturing 
organization where it receives certain Inputs and after 

processing them gives certain outputs for implementation. The 
quality of these outputs determines the effectiveness of a MRP 
system. For MRP to generate proper output, accurate and timely 
availability of the inputs is essential. Accuracy of these Inputs 
greatly affects the performance of MRP. Below we discuss important 
decision making areas in a MRP system. These are lot sizing, 
sequencing, determination of planned lead times and the capacity 
planning. 

1.2 The Lot Sizing decision 

In the lot sizing procedure the timing and size of the 
replenishment order is decided given the discrete period demands 
i.e requirements are converted into a series of replenishment 
orders. For one component this problem involves determination of 
the grouping of the requirements data into a schedule of 
replenishment orders that minimizes the sum of setup costs and 
inventory carrying costs. In a multistage assembly environment 
the decision regarding lot sizing becomes all the more Important 
as lots of the previous machines affect the lot sizing decisions 
of subsequent operations. Thus obtaining optimal lot size 
solutions becomes more difficult. 

Various lot sizing situations arise and optimal lot 
sizing decision is determined for each case by a separate 
procedure. It was argued above that lot sizing rules for a 
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production system where "n" machines are arranged in series is 
likely to be more complicated than the lot sizing rules for single 
machine case. Even for single machine case situation is different 
when demand is deterministic and when it is probabilistic. The 
case of deterministic demand is further classified into two cases 
such as (i) uniform and (ii) dynamic. Simple EOQ formula gives 
optimal lot size decisions for the case of uniform demand. This 
rule produces "robust" lot sizes and hence can be used even there 
is little uncertainty about the demand figures. The dynamic 
programming based on Wagner Whitin rule [13] is popularly used for 
the case of dynamic demand. When the demand is uncertain, the lot 
sizing rules consider the additional factor of service level. 
Hadley and Whitin [7] give a good classification of these rules. 

Since the situation in MRP is a case of multiple 
machines in series and with many products and demand may be 
assumed to be dynamic (which itself is a simplifying assumption), 
the lot sizing decision is not an easy one, atleast 
computationally. Hence simple lot sizing rules such as EOQ and 
Wagner Whitin rules are repeatedly applied even though it is known 
that they will produce sub optimal solutions. Many lot sizing 
rules that are used in MRP contexts are lot for lot, Silver Meal 
heuristic, Wagner Whitin algorithm, Part Period balancing etc. 

1.3 Scheduling decision 

In a multistage problem where there are a large number 
of processing machines, since processing a lot through a work 
centre takes time, the processing of other parts which use the 
same work centre is delayed and they have to wait in a queue. This 
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delays the availability of assemblies at upper levels and leads to 
increased inventory carrying costs. 


Thus 

the sequencing decision 

affects 

the 

due date 

performance as 

well as 

the inventory 

carrying 

costs 

Various 

sequencing rules used 

in MRP context 

are 

SPT 

(shortest 


processing time first), BDD (earliest job due date first), CRR 
(critical ratio rule), ORIB (number of operations remaining in a 
branch) etc. 

1.4 Planned lead time 

Planned lead times are used by the MRP system for the 
purpose of planning. These lead times are used to determine order 
release dates. The lead time of a manufactured item is made up of 
a number of elements. These are - waiting time in queue, running 
(machining, fabrication, assembly etc) time, setup time, waiting 
time (for transportation, inspection etc). 

1 . 5 Capacity Planning 

Capacity requirements planning is the function of 
determining what capacities will be required, by work centre by 
period, in the short to medium range, to meet current production 
goals. The output of MRP system indicates what component items 
will have to be produced and when, and this output can therefore 
be converted into the capacities to produce these items (Orlicky 
[ 8 ]). 

Some of the performance criteria over which these 
decisions are judged are number of late orders, total aggregate 
lateness, total units lateness, total inventory carrying cost, 
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number of stockouts etc. The performance over these parameters is 
affected by the choice of lot sizing rules, scheduling rules, 
planned lead times and capacity planning decisions. It has been 
shown that lot sizing rules and scheduling rules interact with 
each other and that some lot sizing rules go well with certain 
types of sequencing rules (Biggs[2]). Thus the decision maker is 
actually confronted with the decision regarding the combination 
rather than the individual decision of lot sizing or scheduling. 
We discuss these issues in detail in chapter 2 on the literature 
survey. 

1.6 Scope of Present Work 

Most of the earlier works in MRP have been to determine 
the effects of various rules in different decision categories and 
seek to choose the best combination of several rules. Very few 
works are available that help to improve a given solution. The 
work of Sharma and Potty [11] was an exception where they have 
appropriately modified the inventory carrying costs of the 
products used for lot sizing decisions based on the results of 
scheduling decisions. This had resulted in significant improvement 
of results over the best known combination of decision rule for 
lot sizing (which is Wagner Whitin [13]) and sequencing (which is 
the Critical Ratio rule, see Goodwin and Weeks [6]). Later, 
Yoganandan [14] has modified the priorities in the sequencing 
module of MRP based on their earliness/lateness at an assembly 
centre and obtained superior results than the rules given in 
Sharma and Potty [11]. 

We refine the idea given in Yoganandan [14] to develop 
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the concept of coordination index and develop rules which give 
results that are comparable to the ones given by Yoganandan, but 
give significantly better results on the dimension of total units 
lateness. 

In the present work the new approach of achieving 
coordination amongst the child components attempted to reduce the 
range (in terms of time) of child components reaching a assembly 
centre in a multistage assembly environment. For lot sizing and 
scheduling, standard rules which had shown good results in 
earlier work done were chosen (the details are included in the 
chapter on literature survey). By using the coordination range for 
calculating a coordination index which was later used as a 
priority value for lots, a further improvement in the performance 
of the system was achieved which establishes a relation between 
coordination and performance of the system. Two heuristics using 
this approach were developed and one based on considerations of 
lead time was developed. 

1.7 Reasons for using the iterative heuristic approach 

Mathematical model which takes into account the lot 
sizing of parts at different levels of the product structure and 
sequencing rule for scheduling the parts on each work centre is 
quite complex and the solution procedure quite expensive. The 
heuristic approach based simulation is computationally much more 
attractive. 


1.8 Organization of the Thesis 
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In chapter 2 we briefly give the literature survey and 
in chapter 3 we give the prerequisites, tools and terminology of 
MRP. 

In chapter 4 we describe the new heuristics developed by 
us. In chapter 5 we give the essence of the MRP simulator whose 
program listing is given in Appendix A. Basic skeleton of this 
program was proposed by Dr R.R.K.Sharma in which we have embedded 
the various heuristics. In chapter 6 we analyze the results 
obtained and discuss these results. In chapter 7 we present our 
conclusions and discuss the implication for future research. The 
details of twenty test problems with results and " t" test results 


are given in Appendix B. 



CHAPTER 2 


LITERATURE SURVEY 


The impact of product structure on the performance 
of priority dispatching rules was studied by Russell and Taylor 
[10]. Two type of BOMs were studied. It was found that as BOM gets 
taller with more levels, the probability that the job will finish 
after its due date increases. It was also reported that 
performance of shortest processing time (SPT) sequencing rule 
improves as product structure gets taller. 

The effects of product structure and sequencing 
rule on assembly shop performance was studied by Fry, Oliff, 
Minor and Keong [4]. In this paper an effort was made to study 
the effect of product structure on the performance of fourteen 
selected priority rules in a simulated multistage job shop (MRP 
environment). The product structures chosen can be classified into 
three groups - flat, tall and complex. The results established a 
relation between product structure and scheduling rules in the 
performance of MRP system and specific rules are indeed better 
suited to certain product structures. 

The results of the study indicated that there is a 
strong relationship between product structure and sequencing rule 
performance. Taller structured BOMs tended to be more tardy than 
flat BOMs using the same due date allowance, a finding supported 
by Russell and Taylor [10]. But contrary to to Russell and Taylor 
this study indicated that the performance of SPT does not improve 
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as product gets taller. SPT never performed well for flat BOMs for 
any performance measure. Despite being one of the best rules for 
the single stage job shop, results from this study suggest that 
that SPT is not appropriate for a shop that performs assembly 
operations. 


This study indicates that branch critical ratio 
rule (BCR) performs well on the criterion of mean absolute 
lateness (MAL). 


Dispatching in a multistage job shop where machine 
capacities are unbalanced was studied by Fry, Philipoom and 
Markland [5]. In this paper, the performance of selected 
dispatching rules in an unbalanced multistage job shop was 
evaluated to determine whether the performance was consistent for 
jobs which were routed through the bottleneck and for jobs which 
bypassed the bottleneck. The results indicated that sequencing 
rules which are based on due date tend to be the overall better 
performers for both type of jobs. No rule was found to be superior 
in all cases for all measures. Another result was that the results 
for the jobs which pass through the bottleneck were not the same 
for the jobs which bypass the bottleneck. Also the degree of 
capacity imbalance or due date allowance did not cause a 
significant change in the ranking of the sequencing rule 
performance. 

R.R.K.Sharma and V.S. Potty [11] developed multipass 
heuristics for improved stockout performance. In this work two 
heuristics were used to improve the stockout performance. In them 



lot sizes were selectively reduced in steps to get the desired 
improvement. The criterion for choosing the part for revising the 
lot size was based on the difference in its actual and theoretical 
inventory carrying cost. After choosing the part its holding cost 
was changed to revise the lot size. The results indicated that 
increase in the holding cost reduces the lot size which improves 
the stockout performance but excessive decrease in lot sizes may 
increase the number of setups, and adversely affect the due date 
performance (as more capacity may be lost in setups now). To take 
care of this the modification in lot size was limited to two 
times. The study also indicated that critical ratio rule (CRR) 
consistently performed better on a variety of performance 


criteria. 



MRP : TOOLS AND TERMINOLOGY 


The objective of this chapter is to describe the basic 
concepts of MRP, terms used in MRP and the methodology of MRP. The 
objective is to provide a overview. For finer details one can 
refer to orlicky[8]. 

3.1 Prerequisites 

i) Existence of a Master Production Schedule based on forecast of 
product demand is essential for MRP. It states how many items are 
required and when they are required to be produced or purchased. 

ii) Inventory record for all components with inventory status data 
is required as an input for MRP. 

iii) Each item is required to be unambiguously identified by a 
unique component number. A BOM indicating all the components and 
the structure giving all the details regarding processes involved, 
child assemblies, processing and routing sequence is required as a 
input for MRP. 



Outputs 


Figure 3.1 


MRP Inputs 
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3.2 Terminology for MRP 

3.2.1 Pegging 

For many items the demand is created from many sources. 
The source identity gets lost when the total requirements are 

collected for an item’s material requirements plan. For some 
occasions it may be important to know the status of the status of 
the items intended for certain destinations. For instance, if 
production delays are likely to create a shortage of a particular 
item, it would be helpful to know which assemblies, finished 

products and, ultimately, customer orders may be adversely 
affected. This is accomplished by pegging where an item is 

’pegged’ with information to identify its parent in a given time 
bucket, see 0rlicky[8]. 

3.2.2 Netting 

Netting is the process in which net requirements are 

calculated by subtracting quantity available from the gross 
requirements. Gross requirement figure is based on the demand 
forecast for the item. 

3.2.3 Offsetting 

The determination of the appropriate time for release of 
planned orders with the objective of their timely completion to 
meet their demand is called offsetting. The planned order release 
determined by subtracting the lead time from the data of earliest 
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net requirement which it is expected to satisfy, see Orlicky [8]. 

3.2.4 Lot Sizing 

Lot sizing is based on the principle of economy of 
scale. In it by giving due consideration to ordering cost, ( setup 
cost in case of in house production) holding and other costs 
involved lot sizes are produced in order to minimize the overall 
cost involved in the production of an item. There are many methods 
of obtaining economic lot sizes. A list of nine such methods is 
given in orlicky [8]. Here Wagner Whitin [13] is discussed which 
gives optimal solution for the single machine dynamic demand case. 
This is used in the simulation model in the present work. 

3. 2. 4.1 Wagner Whitin algorithm 

This algorithm uses an optimizing procedure based on 
dynamic programming. It evaluates all possible ways of ordering to 
cover net requirements in each period of the planning horizon. Its 
objective is to arrive at the optimum ordering strategy for the 
entire net requirements schedule. The Wagner Whitin algorithm 
minimizes the total cost of setup and carrying inventory, and is 
used as a standard for measuring the relative effectiveness of the 
other discrete lot sizing techniques, see Orlicky [8]. 

3.2.5 Scheduling 

Scheduling in a MRP system involves decision regarding 
priority for processing of components i.e it decides in which 
order components will be processed once they reach the work 
centre. In the simulation model for the present work scheduling 
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has been used for decisions regarding queue formation in front of 
a machine. There are many scheduling rules for job shop 
environment, a detailed list of which is available in Baker [1], 
In the MRP environment the Critical Ratio Rule (CRR) has been 
shown to do well on many performance criteria, see Biggs [2]. 

3. 2. 5.1 Critical Ratio Rule 
Critical Ratio is given by, 

Critical Ratio = Lead time remaining - Work remaining 

Work remaining 

The job having the least critical ratio is given 
the higher priority. This is based on the principle of processing 
that component first which is needed the most in the next period. 
This due date based rule is shown to perform well with respect to 
all the performance criteria (Biggs [2] and Goodwin and Weeks [6]). 
3.2.6 Explosion 

Explosion is the process of translating product 
requirements into component part requirements, taking existing 
inventories and scheduled receipts into account. 

In the explosion process, the demand data from master 
production schedule is taken as the product requirement. Using 

this data and data from BOM file, component requirements are 


calculated. 



HEURISTIC APPROACHES TO MRP 


In the present work three heuristics have been 
developed. Two of the heuristics (which were developed) were based 
on considerations of coordination range and one on lead time. We 
have also tested the performance of the heuristics developed by 
Potty [9] and Yoganandan [14]. 

In the work that was done by Potty, in the 
heuristic that he developed, it was attempted to improve the 
stockout performance of the system by selectively reducing the lot 
sizes. The criteria for selecting the components for lot size 
modification was based on the difference in actual and 
theoretical inventory carrying costs. The components for which the 
difference exceeded ten percent (this limit is set arbitrarily) 
were chosen for lot size modification. The lot sizes were modified 
by increasing the holding costs of the relevant components. This 
approach was reported to have given better stockout performance 
over the original result (the one without the lot size 
modification) . 

In the work which was done by Yoganandan [14], the 
scheduling rule used was modified for better performance. In this 
approach after the scheduling was done, the assemblies arriving 
later/earlier than fifty percent (again this figure of 50% is 
arbitrary) of the subassemblies arriving at an assembly centre 
were identified, and for the subsequent rescheduling their 
criticallity at the previous processing centres was 
increased/decreased. The rescheduling was attempted for a maximum 


of twenty times. 
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of twenty times. 

In the present work the information from Potty [9] 
and Fry et.al [4] was used for choosing the standard lot sizing 
and scheduling rules. It was reported there that due date based 
scheduling heuristics are overall better performers. The first 
simulation run was done using these standard rules. In the first 
run the coordination range ( a term used for difference in time 
between the latest and earliest lot reaching a assembly centre ) 
for various assembly centres was calculated. The mean of the times 
of lots reaching a particular assembly centre was calculated for 
each assembly centre. The priority of the lots at previous 
processing machines was decided by the latest reaching lot of a 
particular assembly. Depending on the time by which the latest 
reaching lot of a component varied from the mean value for a 
particular machine and whether it was more or less than the mean 
value, its priority was modified at the previous processing 
machines. In the finally developed approach only lots of those 
components were modified whose latest reaching lots had reached at 
a time greater than the mean value. 

* Seven parameters were observed as indicators of 

system performance. These were : 

i) Total aggregate lateness - This is the sum of all the times by 
which the lots get delayed. The lower is the value of this 
parameter, the better is system performance. 

ii) Total setup cost - It is sum of all the costs involved in 
setting up the machines for various operations. 

iii) Total units lateness - This quantity is the sum of the 
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product of units and the time period units by which the components 
get delayed. 

iv) Total inventory carrying cost - This is sum of all the 
inventory carrying costs. 

v) Theoretical inventory carrying cost - The theoretical inventory 
carrying cost is given by the lot sizing model which seeks to 
balance the inventory carrying cost and setup cost while ignoring 
the sequencing aspects of the MRP system. 

vi) Average utilization - For the general machines capacity 
utilization is given by the percentage of the total available time 
for processing when the machine was busy. Average utilization is 
the average of all the machine utilizations. 

vii) Average coordination range - Coordination range for a 
particular assembly centre is the difference between the time of 
latest reaching lot at the assembly centre and the earliest 
reaching lot at the assembly centre. The average coordination 
range is average of the coordination ranges of all the assembly 
machines used. 

In the present work, it was attempted to reduce the 
coordination range and see its effect over all the remaining 
parameters which have been described earlier. 

The five heuristics which were tested are as follows. 

4.1 Potty’s Heuristic ( HI ) [9] 

Let R be the set of parts with difference in actual 
and theoretical holding cost exceeding 10 % and S be the set of 
parts already considered twice in the iterative procedure. Let 
K[j] denote the number of times the part j has been considered in 
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the iterative procedure, and BestSol denote the best solution 
arrived at with respect to total units lateness. The steps of the 
heuristic are described below. 

Step 1: Initialization. Set K[j] to zero for all j. Using 
Wagner Whit in for lot sizing and CRR for scheduling simulate the 
material requirements planning system to prepare the set R and 
initialize the following. 

S <- 0, K[j] 4- 0 for all j, BestSol 4 - infinity 

Step 2*. Select part j from set R which has the highest percentage 
difference in theoretical and actual holding costs. Increase the 
holding cost of part j in steps of 0.2 until the lot size changes. 
Set K [ j ] 4- K [ j ] + 1 and if (K[j]=2) then updates S + {j}. If 
the current solution is better than the BestSol w.r.t to total 
units lateness then update BestSol Current solution. 

Step 3: Simulate and prepare the set R afresh and let R <- R - S 
Step 4: If R = 0 then stop else go to step 2. 

4.2 Yoganandan’s heuristic ( H2 ) [14] 

After the scheduling is done, the assemblies 
arriving later/earlier than the fifty percent of the subassemblies 
arriving at an assembly centre are identified, and for the 
subsequent rescheduling their criticallity at the previous 
processing centre is increased or decreased. The rescheduling is 
attempted for a maximum of twenty times. 

4.3 Heuristic developed ( H3 ) 

4.3.1 The Coordination Index approach 

In this approach, based on the coordination range 
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for a assembly centre and the latest reaching lot of a component, 
coordination index for the lots of the component is calculated. 
This value along with the scheduling rule (CRR in the present 
case) is used for deciding the priority of the lots in the 
subsequent simulation runs. 

Let A be the set of components which undergo 
processing. Let C be the set of coordination indexes corresponding 
to each member of A. Let BestSol be the best solution arrived at 
with total units lateness. 

Step 1: Initialization. Set C[j] to one for all j, where j 
denotes a component no. Using Wagner Whitin algorithm for lot 
sizing and CRR for scheduling (i.e for deciding priority) simulate 
a material requirements planning system. Initialize BestSol e- 
inf inity. 

Step 2: First run - In the first simulation run use the scheduling 
rule with coordination index value as one for deciding the 
priority of lots when they join queues in front of general 
machines. 

Step 3: Next run - In subsequent simulation runs set coordination 
index to 

C[j] <r 1 + (latest reaching lot of j - mean time of 

reaching corresponding to assembly machine at next level) 
Set priority value to 

Priority value of j f- (Product of Critical Ratio and 

coordination Index) * (-1) 

where, Critical Ratio is given by, 

Critical Ratio = Lead time remaining - Work remaining 

Work remaining 
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Step 4: Perform step 3 twenty times. 

In the heuristic H3 we kept the lead times fixed 
and did not vary it across simulation runs. 

4.4 Heuristic based on lead time ( H4 ) 

In this, heuristic along with applying coordination 
index for deciding priority as in H3, planned lead times for the 
lots was revised depending upon the waiting time of the lot in the 
queue in front of the general processing machines. The planned 
lead time increase was proportional to the waiting time in the 
queues. 

4.5 Heuristic based on modifying the priority according to 
coordination index of every lot ( H5 ) 

In this heuristic, coordination index corresponding 
to each lot was calculated (unlike H3 where coordination index 
corresponding to each component was calculated). Rest of the 


method was same as H3. 



CHAPTER 5 


IMPLEMENTATION 

In the present work, in order to test the 
performance of the developed heuristic, a simulated job shop with 
assembly and general processing machines was used. A simulator of 
such a shop was already developed by Dr R.R.K.Sharma in PASCAL. 
This program was debugged and linked with the heuristics developed 
in this work. 


5.1 The Simulated Job Shop 

The simulated job shop had assembly machines and 
general processing machines whose number could be changed 
depending upon the requirement of the problem. A problem consisted 
of all the necessary inputs for a MRP system i.e MPS, BOM (product 
structure), processing and setup times for various machines and 
the cost data which included the holding costs, order costs, setup 
costs, processing costs, cost of component etc. 

The BOMs (product structures) chosen were of 
different types as different type of product structures have been 
found to behave differently with a scheduling rule as reported by 
Fry et.al [4]. While designing the problems the time and cost data 
as well as the demand figures for various components was randomly 
generated. 

This data of the the problem was given as a input to 


the system. 

The simulation carried out was a event based 
simulation. The sequence of events that took place in the 
simulation are described in the paragraphs to follow. 



simulation are described in the paragraphs to follow. 


The first operation was to read all the data from 
the relevant files. This was followed by the explosion process 
where the product requirements from the MPS was translated into 
component level requirements. 

Once explosion was over, based on the data from MPS 
and data of costs, lot sizing was done by using Wagner Whit in 
algorithm. This resulted in economic lot sizes to be produced in 
different periods. 

Now based on the data related to product structure, 
those components were chosen which had no child assemblies. 
Depending on which general machine was required for its processing 
first it was put in the queue in front of that assembly machine. 
The position of the lot in the queue was decided by the priority 
value based on Critical Ratio Rule in the first run and by 
combined effect of critical ratio and coordination index in the 
subsequent runs. After all the lots in front of the present lot 
were processed, the present lot was processed. This was followed 
by a trip to a assembly or general machine depending upon the 
requirement. This process continued till the final finished 
was produced. The general approach for lot routing is as shown by 
the diagram on next page. 
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Figure 5. 1 
Simulation Scheme 
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5.2 Test problem development procedure 

The heuristics tested are based on the lot sizing 
decision and the sequencing decision in the MRP systems. Although 
heuristic H4 tries to incorporate the effects of lead time 
decision also, the capacity decision is totally ignored in the 
developed/considered heuristic. 
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However it was found that the developed heuristics 
were totally ineffective when machine capacity utilization 
variance was high. Hence we attempted to control this factor in 
the test problems developed, and tried to generate test problems 
that kept utilizations across machines at uniform level. This 
required lot of time in developing the test problems. Typically it 
took four days to develop a test problem. 



CHAPTER 6 


RESULTS AND DISCUSSION 

On the simulation model developed, twenty test 
problems were tried. While designing the problems it was tried to 
keep a variety of product structures, as it has been reported by 
Fry el.al [4] that different product structures behave differently 
with different scheduling rules. By doing this it was tried to see 
if the developed heuristics were consistent. The various time 
figures ( processing times, setup times etc ) and cost figures 
(holding costs, processing costs, setup costs etc ) were randomly 
generated. The details of test problems and the results obtained 
(including t test results) appear in Appendix B. 


5.1 Comparison of Results 

For comparing the results of the three heuristics, 
test of significance ( t test ) was performed over the data of 
percentage improvement given by each heuristic over the 
combination of Wagner Whitin rule and CRR rule whose results 
appear under the head "original". The three heuristics for which 
data have been considered are HI, H2, H3. Since H4 and H5 did not 
give any improvement over H3 or H2, their data have not been 
included due to lack of time. 

The twenty problems which were tried have been 
included in Appendix B. The analysis of the results of the three 
heuristics on seven parameters is as follows. 

5.2 Discussion 

On the basis of results generated by performing the 
t - tests on percentage improvement figures on various performance 
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criteria, following can be inferred. 

i) Total units lateness - H2 though gives better results than HI, 
the results are not significant at 5% significance level. Results 
given by H3 are significantly better than both HI as well as H2 
over this parameter at 0.5% significance level. In this case, HI 
though gives better results than original in nine out of twenty 
cases, and inferior results only in six cases, the mean 
improvement is negative. This lends credence to the fact that the 
rules are suitable for certain categories of product structures as 
has been found by Fry et.al [4]. 

ii) Total aggregate lateness - The results obtained by H2 are 
though better than HI, they are not significantly better at 5% 
significance level. H3 gives significantly better performance than 
both HI as well as H2 at 0.5% significance level. 

iii) Total inventory carrying cost - Both H2 as well as H3 give 
better results than HI at 0.5% significance level. H3 though gives 
better results than H2, its results are not significantly 
different at 5% significance level. 

iv) Average utilization - H2 gives better results than HI but 
results are not significantly different at 5% significance level. 
H3 gave significantly better performance than both HI as well as 
H2 at 0.05% significance level. 

v) Average coordination range - HI gave significantly better 
performance than H2 at 0.05% significance level. H3 gave better 
performance than both HI as well H2 at 0.05% significance level. 

vi) Theoretical inventory carrying cost - Except for HI where due 
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to lot size change this parameters change, other two do not show 
any modification in this parameter. Therefore, for this parameter 
t test was not applied. 

vii) Setup cost - For this parameter also, except for HI, other 
two heuristics do not show any change. Therefore, this was not 
considered for t test. 

In all the problem instances we had considered, it 
has been uniformly observed that if the machine utilization 
variance was high then it resulted in all the heuristics HI, H2 
and H3 being ineffective in giving improved results. Hence, 
machine capacity utilization had a substantial effect on the 
efficacy of these heuristics. This has also been the observation 
of Potty [9] and Yoganandan [14] who had prepared problem 
instances where all machines had relatively uniform capacity 
utilizations. Similarly we had included sufficient machine 
capacities to keep machine utilizations at a uniform level, the 
heuristics produced the results which have been documented here. 
In future work one has to consider the effect of lopsided machine 
utilizations and derive heuristics that give results in such 
conditions. 

Potty [9] and Yoganandan [14] had kept the lead 
times constant across the lead times constant across all the 
components and assemblies. We had kept it variable and in 
proportion to waiting times before the machines for each component 
and subassembly. This could be an influencing factor that possibly 
has lead to improved performance over the performances obtained in 
Potty [9] and Yoganandan [14]. Future work could be done to refine 


the lead time heuristics even further. 
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Due to lack of time we could not combine the 
heuristics HI and H3 into one procedure where we sequentially 
apply HI and H3 till improvements are possible. When lot size 
modifications cease to improve results, we could initiate 
sequencing oriented changes to seek improvements and so on. 

Yoganandan [14] had shown that sequencing 
modification based heuristic H2 (as reported here) does not 
improve the performance when faced with flat structures. We have 
identified 11 out of 20 problems attempted here where lot sizing 
modification based heuristic H2 (as reported here) does not 
improve the performance. There could be several reasons for this. 

i) The selection of parts qualifying for holding cost modification 
( parts having 10% difference in actual and theoretical holding 
costs) tends to be a bit arbitrary. 

ii) The product structures of each problem instance could be 
different . 

One may conduct simulations controlling these 
factors to determine the strength of their effects. One of the 
findings is that for heuristic HI the increase in costs was 
associated with reduced coordination index for eight out of twenty 
problem instances. It indicates that coordination index is a 
useful measure of shop performance only if assemblies are "in 
time". If all the assemblies are coordinated "late", then it 
ceases to be a useful measure. 

In this research we have presented a refined 
heuristic that seeks to improve the overall shop performance. 
Since our sample size is limited its generalizability is tentative 
even though its performance is impressive. 
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Varying lead times of each lot (heuristic H4) and 
changing priority of each lot (heuristic H5) does not give 
significantly better results, although we did not document these 
results due to lack of time. 

There is scope to do research to find diagnostic 
measures which tell us which heuristic to use (HI, H2 or H3 or any 
other). In this work we have attempted to determine the efficacy 
of two such measures, i.e, (i) Difference in actual and 
theoretical inventory carrying costs and (ii) Coordination range 
defined as the variance/range of arrival times of different 
subassemblies at an assembly centre, with substantial success. 

It has been recognized in literature that certain 
heuristics are more suited to certain product structure and that 
all the decisions in MRP systems (lot sizing, sequencing, lead 
times and planning machine capacities) interact with each other. A 
detailed investigation could throw more light on these issues. 

Significantly superior results obtained by 
heuristic H3 shows that computationally attractive heuristics have 
promise in giving improved performance. 



CHAPTER 7 


CONCLUSIONS 

One of the primary aims of this research was to 
develop a approach by which coordination among the components can 
be achieved, which can improve the system performance. It was 
observed that by reducing the variance in arrival times of 
components reaching a assembly centre (i.e by ensuring better 
coordination), system performance can be improved. The approach 
ultimately followed was of reducing the time range at a assembly 
centre (referred to as coordination range in the present work) 
which was found to give better results than the two works earlier 
done. 

Out of seven parameters considered, total units 
lateness, total aggregate lateness and average utilizations were 
the parameters for which coordination approach gave significantly 
better results. The coordination approach was also found to have 
improved average capacity utilization. 

Scope for future work 

Some related areas where some future work can be 
done are as follows 

i) Use of modification in Lot Sizing to reduce coordination range 
for better system performance. 

ii) Since scheduling rules are sensitive to product structures, 
effect of product structures on coordination range can be studied. 

iii) It was found that by having a better coordination even the 
average capacity utilization gets improved significantly. Study of 
various factors which can improve the capacity utilization can 


also be under taken. 



also be under taken. 


Iv) In many cases the variability in improvement is high, which 
implies that these rules do not perform well on all product 
structures. Hence, one needs to develop a framework which suggests 
a matching rule to a given product structure. This research has 
gathered sufficient evidence to show that such a framework is 
possible. 



34 


REFERENCES 

1* Baker. K. R, "Introduction to Sequencing and Scheduling" , 
John Wiley and Sons, 1974. 

2. Biggs. J.R, "Priority rules for shop floor control in a 
material requirements planning system under various levels of 
capacity", International Journal of Production Research, 1985, 
23(1), 33-36. 

3. Bowker and Lieberman, "Engineering Statistics", Prentice Hall, 
England. 

4. Fry. Timothy. D, Oliff. Michael. D, Minor. Elliot. D, 
Leong. G. Keong, "The effects of product structure and sequencing 
rule on assembly shop performance". International Journal of 
Production Research, 1989, VOL 27, 671-686. 

5. Fry. Timothy. D, Philipoom. Patrick. R, Markland. Robert. E, 
"Dispatching in a multistage job shop where machine capacities are 
unbalanced", International Journal of Production Research, 1988, 
VOL 26, No. 7, 1193-1223. 

6. Goodwin. J.S and Weeks. J.K, "Evaluating scheduling policies 
in a multilevel assembly system", International Journal of 
Production Research, 1986, 24, 247-257. 

7. Hadley and Whitin, "Analysis of Inventory Systems", Prentice 
Hall, India, 1980. 

8. Orlicky. Joseph, "Material Requirements Planning" , Mcgraw 
Hill Book Company, 1975. 

9. Potty. V.S, "Simulation of heuristics, lot sizing and 
sequencing procedures in MRP", M.Tech dissertation, I.I.T, 
Kanpur, 1990. 

10. Russell. R.S and Taylor. B.W, "An evaluation of sequencing 

rules for an assembly shop", Decision Sciences, Vol 16, 196-212, 

1985. 

11. Sharma. R.R.K, Potty. V.S, "Multipass Heuristic for improved 
stockout performance", In Dr V.Raju (Ed): Design automation and 
computer integrated manufacturing, Tata McGraw Hill Publishing 
Company, New Delhi, 1991, 275-282. 

12. Vollman. Thomas. E, "Production Planning and Control Systems", 
Galgotia Publications Pvt Ltd, 1989. 

13. Wagner. H.M and Whitin. T.M, "Dynamic Version of Economic Lot 
Size Model" , Management Science, October, 1958, 89-96. 

14. Yoganandan.S.R, "Coordinated Scheduling in MRP Systems", 
M.Tech dissertation, I.I.T, Kanpur, 1992. 



APPENDIX A 


SIMULATOR DEVLOPED 

The following is the program in PASCAL which was used for 
simulation. 


program mrp( input, output) ; 
const 

max_mc_no 

max_jio_p f _p r o c e s_o n_any_a s s 1 y 
max_assly_no 

max_subassly__no_into__an_assly 
max_pe r i od_no 
infinity 

max_s tatistic s_p t_no 
computing_points__per_period_in_w_w 
max_no_of_lots 
max_queue__length 

inventory_cost_f raction = 0.12; 
no_of_hours_in_a_year = 8640; 
no_of_hours_in__a_j?eriod = 24; 
collect_statistics_from_time = 10.0 * no_of Jhours_in_a_period; 
lead_time=50. 0; 
simulation__time=inf inity; 

type 

range_of__assly_mcs = 16.. 20; 

ptr_to_event = ^event^record; 

event__type = (collect_statistics_f rom_now_on, 
lot_enters_the_shop, 
lot_enters_the_general_mc_queue, 
lot_enters__the__assly_mc_queue, 
process ing_on_general_mc_ < _s tarts, 
processing_on_assly_mc_starts, 
process ing_on__general__mc__gets_over, 
process ing_on_assly_mc_gets_over, 
lot_enters__the_f inished_product_list ) ; 

process_data = record 

mc_no : integer; 

setup_time, 

processing_time_per_unit , 
setup_cost : real; 
end; 

child_assly = record 

asslyjio, 

qty_to_be_contributed : integer; 
end; 

subassly_data = record 

assly__no, 
qty_av : integer; 


= 15; 

= 5; 

= 22 ; 

= max_assly_no; 

= 50; 

= 40000000; 

= 4; 

= 100 ; 

= max__period_no; 
= 125; 
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assly_no, 
qty_av : integer; 
end; 


lot = record 

assly_no, 

original_size, 

qtyJLeft, 

no_of_agm_v i s i t ed „ 
current_mcjno : integer; 

processing_complete : boolean; 

aggregate_lateness_of_lot , 
units_lateness_of_lot , 
actual_inventory_carrying_cost , 
priority_value, 
percentage_completion, 
required_completion_time, 
ac t ua l_comp 1 e t i on_t i me , 
cpt, 

total_wait_time_in_queues_so_far , 
tpt_completed_so_f ar , 
time__of_entry_into__the_shop, 
total_j~equired_processing_time, 
time_of_joining_the_queue, 
time_when_processing__starts : real; 
end; 

basic_data__about_an_assly = 
record 

assly__no, 

assly_no_af ter_processing , 
parentjio , 
no__o f __ch i 1 d__a s s 1 y s , 
no__of_children_for__assembly , 
no_of __processes_needed , 
qty_f or__parent , 
assly_on_mc_jio : integer; 

actual__inventory_carrying__cost , 
theoretical_inventory_carrying_cost , 
total_setup_cost , 
cost : real; 

lot_sizing__to_be_done , 

immediate ly_gets__into_assly, 

immedi ate ly_needs_j)ro cessing : boolean; 

aps 

: array[l . . max_no_of__proces_on_any__assly * 3] 
of process_data; 

array_of_child_assly 

: array[l. . max_subassly__no_into_an_assly] 
of child_assly; 


end; 



lot_data_about__an__assly = 
record 

assly__no, 

lo t_no_which_wi 1 l__enter_the_shop , 
no — of_lots : integer; 

dol : array [ 1 .. max_no_of_lots] of lot; 
end; 

general_machine = record 

tpt, 

total_idle__t ime , 
cu : real; 
end; 

assly__machine = record 

f inished_assly_no : integer; 

process ing_time_per_unit, 
tpt, 

total_idle_time, 
cu : real; 
end; 

event_record = record 

event : event_type; 

array_of_relevant_nos : array [1.. 3] of integer; 
(* 1 - asslyjno 

2 - lot_no 

3 - mc__no_to__visit *) 

execut longtime : real; 

next_event : ptr_to_event ; 
end; 

(* BEGIN TYPES FOR EXPLOSION DETAILS *) 

explosion_details = record 

assly_no, 

qty_to_be__exploded : integer; 
end; 

(* END TYPES FOR EXPLOSION DETAILS *) 

(* BEGIN TYPES FOR QUEUES *) 

lot__and_assly__no = record 

lot_no, 

assly__no : integer; 
end; 

type_f or_queue_before__general_mc 

= record 

queue_length : integer; 

detail s_of_queue : array [1 . . max__queue___length] 

of 1 o t _and_as s 1 y_no ; 



end; 
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type_for_queue_before_assly_mc = 
record 

f ini sh_ass ly_no , 

no_of_subasslys_into_assly : integer; 

quant ity_ava liable : array [ 1 . . max_subassly_no__into__an_assly] 

of integer; 


details_of_queues 

end; 


array [ 1 . . max_subassly_no__into__an_assly] 
of type_f or_queue_bef ore_general__mc; 


(* END TYPES FOR QUEUES *) 


(* BEGIN TYPES FOR WAGNER WHITIN ALGORITHM *) 

type basicjperiod_data = record 

order_cost , 
holding_cost : real; 
demand : integer; 
end; 

input_data_for_lot_sizing = array [ 1 .. max_period_no] of basic__period_data; 

orderjdata = record 

period__no, 

qty_produced : integer; 
end; 


computation_data = record 

cost: real; 

order jperiod_no : integer; 
end; 


output_data__of_lot_sizing = 
record 

no_of_orders_placed : integer; 

order_schedule : array [ 1 .. max_period_no] of order_data; 
total__cost , 
holding__cost , 
order_cost : real; 
end; 

period_wise_computations = 
record 

no^of^calculations, 
min_cost_position : integer; 

period_wise_calculation : 

array [1. . computing joints_per_period_in_w_w] of 
computa t i onjdat a ; 

end; 

type_jfor__f ind_min_position = array [ 1. . max_period_no] of real; 

(* END TYPES FOR WAGNER WHITIN ALGORITHM *) 
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var clock_time :real; 

demand , 

exploded__demand : arrayfl. . max_assly__no , 1 . . max__period__no] of Integer 
event_freetop, 

ptr_to__current__executable__event „ 
top_of_event__list , 
new___event , 

P : ptr__to_event ; 

abda 

: array [1 .. max__assly_no] of basic_data_about__an_assly ; 

ala 

: arrayfl. ,max__assly_no] of lot_data__about_an_assly; 
agm : 

array [1. . max_mc_no] of general_machine; 

aam : 

array [range__of__assly_mcs] of assly__machine; 

list__of_end_products : array [ 1 .. max_assly_no] of integer; 

(*BEGIN VARIABLES FOR EXPLOSION DETAILS *) 
no_of_asslys_on_explosion_list , 
no__of__end_products : integer; 

incremental_explosion_list : 

array [1. . (50*max_assly_no ) ] of explosion__details; 

(* END VARIABLES FOR EXPLOSION DETAILS *} 

(* BEGIN VARIABLES FOR QUEUES *) 

queue_bef ore_general_mc : array [1. . max__mc__no ] of 

tyP e _f° r _queue_bef ore_general_mc; 

queue — bef ore__assly_mc : array [range_of_assly_mcs] of 

typ e_f or__queue_bef ore_asslyjmc ; 


(* END VARIABLES FOR QUEUES *) 

(* BEGIN VARIABLES FOR WAGNER WHITIN *) 

input_dataJfor_w__w : input_data_f or_lot__sizing; 

output_data__of__w_w : output__data_of_lot__sizing; 

computation__table__forjtf_w: array [1. . max_jperiod_no] of 

period_wise__computations ; 

array_for__output__of_w_w : array [1 .. max_assly_no] of 

output_data__of_JLot_sizing; 

in_ass_dat : text ; 
in_dem_dat : text ; 
in_ass_mc : text ; 
w_w_output : text ; 
w_w_file : text; 
event_JTile :text; 
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potty_output : text ; 
s__dat : text; 
i,j : integer; 
counter: integer; 
in_out :text; 
rand :real; 
flag : boolean; 

mcjDUsy : array [ 1. .maxjno_of jroce.s_on_any_assly] of boolean; 

(* END VARIABLES FOR WAGNER WHIT IN *) 

(* start variables for potty's algo *) 

R,K,S: array [1 . . max_assly_no] of integer; 

theoretical__cost , actual_cost: array [1 . . max__assly_no] of real; 

best_sol : real ; 

count , max_count : integer ; 

stop, go_ahead: boolean; 

total__units_lateness: real ; 

temp, max__dif : real ; 

index, period: integer; 

old_prders,new__orders: integer; 

assly__af ter_jprocessing: integer; 

iteration: integer; 

levels: integer; 

coordination_index, 

present_coordination_index, 

ag_coordination__index : array [1. . max_assly_no] of real; 

max_t ime , max_index , 

present__mean__time, 

prevjnean_time : array [range_of_asslyjncs] of real; 

critical_child : array [ 1 .. max_asslyjno] of integer; 
times: integer; 

arr ival_t ime: array [1 . . 100, 1. . max__assly_no , 1. . max_no_of_lots] of real; 
check_out : text ; 
a,b, c, lots: integer; 

assly_count : array [range__of__assly_mcs] of integer; 

time_of_reaching :array[l. . max_assly_no, 1 . . max_no_of_lots] of real; 

latest_assly__time : array [1. . max_assly_no, range_of_assly_mcs] of real; 
f irst_run : boolean; 
earl iest_lot_t ime, 

latest_lot_time : array [range_of_assly_mcs] of real ; 
prev_range : real ; 


procedure event_addf reelist (p : ptr_to_event ) ; 

begin (* event_addf reelist *) 

p~ . next__event : =event_f reetop; 
event_freetop: =p; 

end; (* event_addf reelist *) 

procedure newmodif ied(var p : ptr_to_event ) ; 

begin (* newmodif ied *) 

if (event_f reetop = nil) 
then new(p) 
else begin 

p: =event_f reetop; 

event__f reetop: =event__f reetop"' . next_event ; 



end; 

p~ . next_event : =ni 1 ; 
p" . array_of_relevant__nos [ 1 ] : =0 ; 
p~ . array_of_relevant__nos[2] : =0; 
p" . execution_time: =0.0; 
end; (* newmodified *) 


procedure print_the_event_list; (* this is the changed version 
var 

q : ptr_to_event ; 


begin (* print_the_event_list *) 

q: =ptr_to_current_executable_event ; 
case q~ . event of 

collect_statistics_from_now_on : 

write ( even t_Jf ile, * collect_statistics_f rom_now_on * ) ; 
lot_enters_the_shop : 

write (event_f ile, ’ lot_enters_the_shop ’ ); 
lot_enters__the_general_mc_queue : 

write (event_f ile, * lot_enters__the_general_jnc_queue 9 ) ; 
lot_enters_the_assly_mc_queue : 

write (event__f ile , 9 lot_enters_the__assly__mc_queue ’ ) ; 
process ing_on_general_mc_s tarts : 

write (event_f ile, ’ processing_on_general_mc_starts * ) ; 
process ing__on__assly__mc__s tarts : 

write (eventjf ile, * process ing_on_assly__mc_s tarts ’ ) ; 
process ing_on__general_mc_gets_pver : 

write (event_f ile, * processing_on_general_mc__gets__over * ) ; 
process ing_on_assly_mc_gets_over : 

write (event_f ile, ’ processing__on_assly_mc__gets__over * ) ; 
lot_enters_the__f inished_product_list : 

write (event_f ile, * lot__enters_the_f inished_product_list’ ) 
end; (* case *) 

write (event__f ile, * assly_no: ’ ,q" . array_of_relevant_nos [1 ] : 1 , * 

’ lot_no: * , q" . array_of_relevant jnos [2] : 1 , ’ * ) 
writeln(event_f ile, ’ exec_time’ , q" . execution_time: 8: 3) ; 
end; (* print_the_event_list *) 


procedure open_j?rint_the_event_list ; 
begin 

rewrite (event_f ile , 9 event, file’ ) ; 

end; 

procedure close_pr int_the_event_l 1st ; 
begin 

close (event_f ile) ; 

end; 

procedure lot_sizing_by_w__w(assly_no : integer); 
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var periodjio, 

order_can_be__placed_f rom_period_no : integer; 

procedure f orward_recursion_of_w_w; 

var period__no, 
index, 
min__pos, 
scratch, 
size, 

no__of__calculations_done, 
search_index_in__previous_period, 
min_cost_pos_upto__last__period, 
order jplaced_in_period : integer; 

min_cost_upto_last_period, 

cost, 

min_cost_upto jeriodJLn_which s _order_is_placed : real ; 

a : type_f orjf ind__min_joosition; 

procedure f ind_min_position(a: type_f orjf ind__minjposition; 

size: integer; 

var minjpos : integer); 


var index : integer; 

min_cost : real; 

begin (* f ind__min_position *) 
min_cost:=a[l] ; 
min_pos: =1 ; 

for index: =2 to size do 
if (a [index] < min_cost) 
then begin 

min_cost : =a [ index] ; 
min_j)os: =index; 
end; 

end; (* f ind_minjposition *) 

begin (* f orward_recursion_of_w_w *) 

computation_table_for_w_w[l ] . no_of_calculations: =1 ; 
computation_table_f or_w_w[l ] . min_cost__position: =1 ; 
computat ion__table__f or_w_w [ 1 ] . per iod_wise__calculat ion [ 1 ] . cost : = 
input_data__for_w_w[l ] . order_cost; 
computation_table_for_w_w[l] .period_wise_calculation[l] . 
order_per iod_no : =1 ; 

order_can__be_jplaced_f rom__period_no : =1 ; 
for periodjno: =2 to maxjperiod_no do 
begin 

no_of__calculations_done: =1 ; 
computation_table__for_w_w [period_no] . 

no_of_calculations:=period_no - order_can_be_placed_f rom__period_no+l ; 
order_placed_in_period: =period_no ; 
min_cost_pos_upto__lastjperiod: = 
computat ion_table_f or_w — w [period_no-l ] . min_cost_position; 
min_cost_upto_last_period: = 
computat ion_table_f or_w_w[period_no-l ] . 
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period_wise_calculation[min_cost__pos__upto__lastjperiod] . cost; 
computat ion — table_Jfor_w_w[period_jio] . 
period__wise_calculation[l ] . cost : = 
min_cost__upto_last_period + 
input __d at a_f o r_w_w [ p e r i o d_n o 3 . order_cost; 
computat ion_table__for__w_w [per iodjno] . 
period_wise_calculation[no__of_calculations_done] . 
order_period_no: =orderjplaced_in_jDeriod; 
search__index_in_previous_period: =1 ; 
order_placed_in_j>eriod:=order_placed_in_period - 1; 

while ( order_placed_in__per iod >= order_can_be_placed_f rom period no ) do 
begin 

no __ 0 f- ca l cu l at i°ns_done: =no_of_calculations_done + 1; 
min_cost_upto_jDeriod__in_which__order__is_p laced: = 
computation_table_for_w__w[period__no--l ] . 

period_wise_calculat ion [search__index_injprevious_period] . cost ; 

cos t : =: min cost up to period_in_jwhich order is placed ; 

for index := order jplaced_in__per iod to (per iod_no-l ) do 
cost:=cost + input_data_for__w_w [index] .ho lding__cost * 
input__data_for_w_w[period_jrio] . demand; 
computat ion_table_jfor_w_w [ per iod_no] . 

period_wise__calculation[no_of_calculations_done] . cost : =cost ; 
computat ion_table_f or_w_w[period_no] . 

period_wise_calculation[no__of__calculations_done] . 
order_period_no: =order__placed_in_j)eriod; 
search_index__in_previous__period: =search_index_in_previous_period + 1 ; 
order_placed__in_period: =order_placed_in_period - 1; 
end; (* while *) 

size:= computat ion_table_f or_w_w [period_no ] . 

no__of_calculations ; 
for index: =1 to size do 

a [index] : =computation_table_for_w_w[period_no] . 

period___wise__calculation[ index] . cost; 
f ind_min_position(a, size, min__pos) ; 

computation_table___for_w_w[period_no] . min_cost_position: =min_pos; 
order_can__be_placed_f rom_period_no: = 
computation_table_for_w_w [period_no] . 
period_wise_calculation[min__pos] . 
end; (* for *) 

end; (* f orward_recursion_of_w_w *) 

procedure f i 1 l_up_the__order_schedule ; 

var prodn_f rom_period, 

period_no, 
min_position, 
no_of_orders_placed, 
index, 

order_qty : integer; 

order_cost , 
total_cost , 

holding_cost : real; 

computat ions_to_continue : boolean; 
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begin (* f ill_up_the_order_schedule *) 

min__position: =computation_table_for_w_w[max_period__no] . 

min_cost_position; 

total_cost : =computation_table_for_w_w[max__period_no] . 

period_wise_calculation[min__position] . cost; 

order_cost : =0 . 0 ; 
no_of_orders_placed: =0; 
period__no: =max_jperiod__no ; 
computat ions_to_cont inue : =t rue ; 
while ( comput at ions_to_cont inue ) do 
begin 

min_position: =computation_table_for_w_w[period_no] . 

min_cost__position; 

prodn_from_jperiod: =computation_table_for_w_w [period_no] . 

period__wise___calculation [min_position] . 
order_period_no ; 

order_cost : =order_cost+input__data_f or_w_w [prodn_f romjperiod] . order_cost ; 
no_of_orders ^placed: =no_of_ordersjplaced + 1; 
order__qty : =0; 

for index: =prodn_f romjper iod to period___no do 

order_qty : =order__qty + input_data_for_w_w[ index] .demand; 
output__data_of__w__w. order_schedule [no__of_orders_placed] . period_no : = 
prodn_f rom_jperiod; 

output__data__of_w_w. order_schedule [no_of_orders_p laced] - qty_produced: = 
order_qty; 

period_no: =prodn_from_period - 1; 

if (period__no = 0) then computations__to_continue: =false; 
end; 

holding_cost : =total_cost - order__cost ; 

output^ data_of_w_w. no_of_orders_placed: =no_of_orders_p laced; 
output__data_of__w__w. total_cost : =total_cost ; 
output__data__of__w__w. holding_cost : =holding_cost ; 
output_data_of__w_w. order_cost : =order_cost; 

end; (* f ill_up_the_order_schedule *) 

procedure print_w__w__output ; 

var index, 

period_no, 

n o _ o f_items__printed_s o _fa r , 
no_of_i tems_to_be _jpr inted , 
no_of_calculations_done : integer; 


begin (* print_w_w_output *) 

writeln(w__w__output assembly no.= ’ , assly_no: 2) ; 

for period_no:=l to max_period_no do 

begin 

writeln(wjw_output , * period_no : = * , period_no: 2, * * ); 
no__of_calculations__done : =computation_ table_for__w__w[period_no] . 

no^of^calculat ions ; 

no_of_items__to__be__pr inted: =no__of_calculations_done; 
no_of_items_printed_so_far : =0; 
while (no_of_items_to_be_printed > 0) do 
if (no___of_items_to_be_jpr inted <= 5) 
then 
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begin 

write (w_w__output , ’ ’ : 12) ; 

for index: = (no_of_items_printed__so_far+l) 

to 

no__of__calculations_done 

do 

begin 

write (w_w__output , computation_table_f o r _w_w [ pe r i o d__no ] . 

period_wise_calculation[ index] . cost:8:3, ’ (’ ) ; 
write (w_w_output , computation__table__for__w_w[period_no] . 

period__wise_calculation[ index] . 
or derjeriod_no, * ) * ) ; 

no_of_items_to_be_printed: =no_of_items_tojDe_printed-*l ; 
end; (* for index *) 
writeln(w_w__output ) ; 

no_of_items__printed_so_far : =no_jDf_calculations_done; 
end (* then *) 
else 
begin 

write (w_w_output , ’ * : 12) ; 

for index: = (no__of_items__printed_so_far+l ) 

to 

(no_of_items_printed_so_far+5) 

do 

begin 

write (w__w__output , computation_table_for_w__w[period_no] . 

period_wise_calculation[ index] . cost : 8: 3, ’ ( * ) ; 
write (w_w_output , computation_table_for_w_w[period__no] . 

period_wise_calculation[ index] . 
order jperiod__no, ’ ) ’ ) ; 

end; 

writeln(w_w__output ) ; 

no_of_items_j)rinted_so_far : =no_of_items_printed_so_far+10; 
no_of_items_t o _be_printed: =no_of_items_to_bejprinted-10; 
end; (* else *) 
end; (* for period_no *) 
writeln(w_w_output ) ; 
write In (w_w__output ) ; 
writeln(w_w__output ) ; 

writeln(w__w_output product ion__schedule is as follows’ ); 
writeln(w_w_output , * period^no’ , ’ production_quantity’ ); 
for index :=1 to output j3ata_of_w_w.no_of_prders_placed do 
writeln(w_w_output , output__data_of_w_w. order_schedule [ index] . 

period_no: 8, ’ ’ : 5, 

output_data_of__w_w. order_schedule [ index] . 
qty_produced) ; 

writeln(w_w__output , ’ order_cost : =’ , output_data_of_w_w . order_cost ) ; 
write In (w_w_output , ’ holding_cost : =* , output_data_of_w_w . holding_cost ) ; 
writeln(w_w_jDutput , 9 total__cost: =’ , output__data_of_w_w. total__cost) ; 
end; (* print_w_w__output *) 

begin (* lot_sizing_by_w_w *) 

forward_recursion_of_w_w; 

f i 1 l_up___the_order_schedule ; 

( * pr int_w_w__output ; * ) 
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end; (* lot_sizing_by_w_w *) 


procedure read_assly_data; 

var 


assembly_no, 
child__assly_no , 
needs jProcess , 
needs__assembly, 
process_no : integer; 

begin (* read__assly_data *) 

reset (in_ass_dat , * in_ass.dat’ ); 

for assembly_no: =1 to max_ assly_no do 

begin 

readln( in_ass_dat , 

abda[assembly_no] .assly_no, 
abda[assembly_no] . 

no_of_chi ld_asslys , needs__process , needs_assembly , 
abda[assembly_no] . cost) ; 
abda[assemblyjao] . 

immediately_needs_processing: = (needs_process = 1); 
abda[assembly_no] . 

immediately_gets_into__assly : = (needs_assembly=l ) ; 
if (abda[assembly_no] . 

no__of___child_asslys > 0) 

then 

for child__assly_no: =1 to 
abda[assembly_no] . 

no_of_child_asslys do 

readln ( in_ass___dat , 
abda [assembly__no] . 

ar ray__of __chi Id _ass ly [ chi Id _ass ly_no ] . ass ly_no , 
abdat assembly jno] . 

array__of_child_assly [child_assly_no] . qty_to_be_contributed) ; 
if (abdatassemblyjno] . immediately_needs_processing) 
then 
begin 

readln ( in _ass _dat , 
abdatassemblyjno] . 

no_of_childrenJfor_assembly, 
abdatassemblyjno] . 

asslyjno_af terjDrocessing, 

abdatassemblyjno] . 

n°jOf jProcesseSjneeded) ; 

for process jno: =1 to 
abdatassemblyjno] . 

no jOf jProcesseSjneeded do 

begin (* for *) 

with abdatassemblyjno]. 

aps [process jno] do 
begin (* with *) 
readln ( in_asSjdat , 

mc_no, setup jtime, processing jtimejPerjUnit, 
setup jCOst ) ; 
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end; (* with *) 
end; (* for *) 
end 

else (* collect_parent_data *) 
begin 

abda[assembly__no] . 
n°_of jprocesses_needed : =0 ; 
readln ( in_ass_dat , 

abda[assembly_no] . parent_no, 
abda[assembly_no] . 

qty_forjparent , 

abdalassemblyjno] . 

assly_on__mc_no) ; 

end; 

end; (* for *) 
end; (* read_assly__data *) 
procedure read_assly_demands; 
var 


assly_no, 
periodjno, 
no_of _demands , 
demand_no : integer; 

begin (* read_assly_demands *) 

for assly_no:=l to max_assly_no do 
for periodjno: =1 to maxjperiodjio do 
demand [ ass lyjio, per iodjio] :=0; 
reset ( in_dem__dat , ’ in_dem.dat* ) ; 
while not eof ( in_dem_dat ) do 
begin 

readln ( in__dem_dat , asslyjio , no jDf ^demands ) ; 
for demand_no: =1 to no_of ^demands do 
begin 

read ( in__dem_dat , per iodjno ) ; 

readln (in_demjdat, demand [asslyjio, per iodjio] ) ; 
end; 
end; 

end; (* read_assly_demands *) 

procedure read_input_data ; 

begin (* read_input_data *) 

read_assly_data; 
read_asslyjdemands; 

•end; (* read_input_data *) 

procedure transf er_asslysj)n_explosion_list (periodjno : 

var asslyjio : integer; 


integer) ; 
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begin (* transfer_asslys_on_explosion_list *} 
no__of __as s 1 y s_on_exp 1 o s i on_l i s t : =0 ; 
for assly___no: =1 to max_assly_no do 
if (demand [assly_no, per iod_no] > 0) 
then 

begin 

no_o f _a s s 1 y s__o n_exp losion__list: = 
no__of_asslys_on_explosion_list + 1; 
incremental__explosion_list [no_of__asslys_on__explosion_list ] . 
assly_no: =assly_no; 

incremental_explosion_list [no_of_asslys_on_explosion_list ] . 

qty_J:o__be_exploded : = 
demand [ ass ly_no, per iod_no] ; 
end; 

end; (* transfer_asslys_on_explosion_list *) 

procedure do_explosion; 

var current_no, 

current_assly_no , 

current_qty_f or__explosion, 

no_of_chi ldren, 

children__no, 

ending_no, 

period_jno, 

temp_assly, 

temp_qty, 

a,b : integer; 

begin (* do__explosion *) 

for a:= 1 to max_assly_no do 

for b:= 1 to max__period_no do 

exp 1 oded_demand [ a , b ] : =0 ; 

for period_no:=l to maxj>eriod_no do 

begin 

transfer_asslys_on_explosionJList (period_no) ; 
current_no: =1 ; 

ending_no: =no_of_asslys_on__explosion_list ; 
while (current__no <= ending_no) do 
begin 

current__assly_no : =incremental__explosion_list [current_no] 

. assly_no; 

current__qty__f or_explosion: = 

incremental_explosion_list [current_no] . 

qty_to_be_exploded; 

exploded_demand[current_assly_no,period_no] : = 

exploded_demand[current_assly_no, period_no] 
+current_qty_for_explosion; 

no_of_chi ldren: = 

abda [current_assly_no] . 

no_of_chi ld_asslys ; 

for chi ldren_no :=1 to no_of_chi ldren do 
begin 

no_of_asslys_on_explosion_list : = 

no_of_asslys_on_explosionJList + 1; 
temp_assly : =abda [current_assly__no ] . 
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array_of __chi ld__assly [ chi ldrenjio ] . assly__no; 
temp_qty: =abda[current_assly_no] . 

array__of_chi ld__assly [ chi ldren_no ] . qty_to__be__contr ibuted 

* current_qtyjEor_explosion; 

exp 1 oded_demand [ temp_assly , per iod_no ] : = 

exploded_demand [ temp_assly , period^ no ] + current_qty__f or__explosion ; 
end ; 

current jno : =current_no + 1 ; 
end; 
end; 

end; (* do_explosion *) 


procedure do__l o t_s i z i ng ; 

var setup_cost : real; 
assly__no, 

assly_no_af ter__processing, 
period_no, 

processjno : integer; 


begin (* do_lot__sizing *) 

rewrite (w_w_output , ’ w_w. out’ ) ; 
for asslyjno : =1 to max_assly_no do 
abda[assly__no] . 

1 o t _s i z i ng_t o__be_done : = t r ue ; 
for assly_no:=l to max__assly__no do 
begin (* for *) 

if (abda [asslyjno] . 

immediate ly_needs_processing) 

then 

begin (* then *) 

assly_no_af ter ^processing: = 
abda [ass lyjio] . 

ass ly_no_af ter__processing ; 
abda[assly_no_af ter__processing] . 

lot_sizing_to_be_done: =false; 
end; (* then *) 
end; (* for *) 

for assly_.no: =1 to max_assly_no do 
if (abda [assly_no] . 

lot_sizing_to_be_done ) 

then 

begin (* then *) 

for period__no: =1 to max_period_no do 
input_data_forjtf_w[period_no] . holding_cost: = 
abda [asslyjno] .cost* 
inventory_cost jfraction; 
setup_cost : =0 . 0 ; 
for processjio: =1 to 

abda [asslyjno] . 

no_of_processes__needed do 

setup_cost : =setup_cost + 



abda[assly_no] . 

aps[process_no] . setup_cost; 
abda [assly__no] . total_setup_cost : = 
setup_cost ; 

for per iod__no : =1 to max_jperiod_no do 
input__data_f or_w___w [per iod__no ] , order__cost : = 
setup^cost; 

for period__no: =1 to max__period_no do 
input_data_for_w_w[period_no] . demand: = 
exploded_demand [asslyjno , per iod_no ] ; 
lot__sizing_by__w_w(assly_no) ; 

array_for_output_of_w_w[assly_no] : =output_data_of_w_w; 
if ( abda [ ass ly__no] . 

immediately__needs_processing) 

then 

begin (* then *) 

assly_no_af terjprocessing: = 
abda [asslyjno] . 

assly_no__af ter_processing; 

array JFor_output_of_w_w[assly_no_after_processing] : = 
output_data_of_w_w; 
end; (* then *) 
end; (* then *) 
close (w_w_output ) ; 

end; (* do_lot_sizing *) 


procedure prepare^ queue_of__lots ; 

var (* lead_time, *) 
assly_no, 
no_°f_lots, 
processjao, 

lot_no : integer; 

begin (* prepare__queue_of__lots *) 

for asslyjno :=1 to max_assly_no do 
begin 

ala [asslyjno] . assly_no: =assly_no; 
ala [asslyjno] . no__of_lots: = 

array_forjDutputj3f_w_w[asslyjno] . no_of_orders_placed; 
no_of_lots : =array_for_output_of_w__w[assly_no] . 

no_pf_orders_placed; 
for lot__no:=l to no_of_lots do 
begin 

ala [asslyjno] . 

dol [lotjno] . asslyjno: =assly_no; 
ala[asslyjno] . 

dol [lotjno] . original_size:= 
ar ray_f or jnutput jnf _w_w [ asslyjno ] . 

order jschedule [ lotjno ] . qty_produced; 
ala [asslyjno] . 

dol [lot_.no] . qty_left: = 
ar ray Jf or jnutput jnf jjw [ asslyjno ] . 

order_schedule [ lot jno ] . qty_produced ; 



ala [assly_.no] . 
dol [ lot_.no] . 
no_of_agm_visi ted: =0; 
ala[assly_no] . 
dol [ lot_no] . 
current_mc_jio : =0 ; 
ala [assly_.no ] . 
dol [lot_no] . 

processing_complete: =false; 
ala[assly_no] . 
dol [lot_no] . 

actual__inventory_carrying_cost : = 0 . 0 ; 
ala[assly_no] . 
dol [lot_no] . 

priority_value: =0. 0; 
ala [assly_.no ] . 
dol [ lot_no] . 

percentage_completion: =0.0; 
ala[assly_no] . 
dol [ lot_.no] . 

required_completion_time : = 
array_for__output_of_w_w[assly_no] . 

order_schedule [ lot__no] . period_no* 
no_of_h°urs_.in_a_period; 

ala[assly_no] . 
dol [ lot_no] . 

actual_completion_time: =0 .0; 
ala [assly_no] . 
dol [lotjno] . 
cpt : =0.0; 
ala[assly_no] . 
dol [ lot_no] . 

total_waIt_.time_in_.queues__so_far : =0.0; 
ala[assly_no] . 
dol [ lot_.no ] . 

tpt_completed_so_f ar : =0 . 0 ; 
ala [assly_.no ] . 
dol [ lot_no] . 

t ime_pf_entry_into_the_shop : = 
ar ray_f or_output_of _w_w [ ass ly_no ] . 

order_schedule [ lot_no] . period_no*l . 0 - lead_tlme; 
ala [assly__no] . 
dol [ lot_no] . 

total_required_processing_time: =0. 0; 

for process_.no: =1 to abda[assly_no] . 

no_of__processes_needed do 

begin 

ala [assly_.no] . 
dol [ lot_.no] . 

total_required__processing_time: = 
ala [assly_.no] . 
dol [lot_no] . 

total._required_processing_.time + ^ 

abda[assly_no] . \ 

aps [process_.no ] . 

setup_time + 4^ Ma. A 

abda [ assly_.no ] . 
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aps [processjno ] . 

process ing_time_per_unit * 
array_for_output_of_w_w[assly_no] . 
order_schedule [lot__no] . qty_produced; 

end; 

ala[assly_no] . 
dol [lot_no] . 

time_of_joining_the_queue: =0. 0; 

end; 

end; 

end; (* prepare__queue_of_lots *) 


procedure plug_in_the_event (p : ptr_to__event ) ; 

var execution_time__just_greater , 

prev_to_execution_time_just_greater : ptr_to_event ; 

begin (* plug_in_the_event *) 
counter : =counter + 1 ; 
if (top_of_event_list = nil) 
then 

begin (* top__of_event_list = nil *) 
top__of_event_list : = p; 
p: =nil ; 

end (* top_of_event_list = nil *) 
else 

begin (* contains at least one event *) 

execut ion_time_Just_.gr eater: =top__of_event_list ; 
prev_to__execution__time__just__greater : =top_of__event__list ; 

if ( top_of_event_list~ . execut ion__time < p~. execut ion__time) 
then 

begin (* continue search till execut ion__time_just_greater *) 
execut ion_t ime_Just_greater : = 

execut ion_time_just_greater~ . next_event ; 

while ( execut ion_time_just_greater <> nil) 

and 

(execution_time_just_greater^ . 

execut ion__time < p" . execution_time) 
do 

begin (* while *) 

execution__time__just__greater : = 

execut ion_time_just_greater* . next_event ; 
prev_to_execution_time_just_greater : = 

prev_to_execution_time_just_greater^ . next_event ; 
end; (* while *) 

if (execution_time_Just_greater <> nil) 
then 

begin 

p" . next_event : = 

prev_to_execution_time__just__greater~ . next__event ; 
prev_to__execution_time_Just_greater" . next__event : =p; 

p: =nil ; 
end 
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else 

begin 

prev_to_execution_time_just_greater'' . next event: =p; 

p: =nil ; - F ’ 

end; 

end (* continue search till execution^ ime_just_greater *) 
else 

begin (* execution time just greater *) 
p~ • next_event : =top_of_event_list ; 
top_of_event_list: =p; 
p:=nil; 

end; (* execution time just greater *) 
end; (* contains at least one event *) 

end; (* plug_in_the_event *) 

procedure prepare_events_regarding_entry_of_components ; 

var assly_.no , 

no_of_lots , 

lot_no : integer; 

begin (* prepare_events_regarding_entry_of_components *) 

for assly_no:=l to max_assly_no do 
if (abda [assly_.no ] . 

no_of_child_asslys = 0) 

then 

begin (* then it is a component assly *) 
no_of_lots: =ala [assly_no ] . no_of_lots; 
for lot_.no: =1 to no_of_lots do 
begin (* for lot_no * ) 
newmodif ied(p) ; 
with p~ do 
begin 

event : =lot_enters_the_shop ; 
array_of_relevant_nos [ 1 ] : =assly_no ; 
array_of_relevant_nos [2] : =lot_no ; 
execution_time: =ala [assly_no] . 

dol [ lot_no] . 

time_of_entry_into_the_shop; 

plug_in_the_event (p) ; 
end; 

end; (* for lot_no * ) 
end; (* then it is a component assly *) 

end; (* prepare_events_regarding_entry_of_components *) 
procedure prepare_collect_statistic_event; 
begin (* prepare_collect_statistic_event *) 

newmodif ied(p) ; 
with p~ do 
begin 

event : =col lect_statistics_f rom_now_on; 
array_of_relevant_nos [ 1 ] : =0 ; 
array_of_relevant_nos [2] : =0; 



execution_time : =collect__statistics_from_time; 

plug_in_the__event (p) ; 
end; 

end; (* prepare_col lect_statist ic_event *) 

procedure initialisations ; 

var mc_no , 

assly_mc_no, 

assly_no_af terjprocessing, 

f inished_assly_no, 

asslyjno, 

subasslyjno, 

quelg_no : integer; 

in_ass_mc : text; 

begin (* initialisations *) 

clock_time: =0.0; 
counter : =0; 
flag: =f alse; 

(* now all me queues and assly me queues are initialised *) 

for mc_no:=l to max_mc_no do 
begin 

queue__bef ore_general__mc [mc_no] . queue_length: =0; 
for quelg__no: =1 to max_queue_length do 
begin 

queue_before_general_mc [mcjio] . 

detail s_o f _queue [ que 1 g_no ] . lot_no:=0; 
queue_bef ore_general__mc [ mc_no ] . 

details_of_queue [quelg_no] . assly_no: =0; 

end; 

end; 

reset (in_ass__mc > * in_ass.mc’ ) ; 
for assly_mcjno: =16 to 20 do 
begin 

readln ( in__ass_mc > 

aam [assly_mc_no ] . 

f inished_assly_no , 

aam [assly__mc_no ] . 

processing__time_j)er__unit , 
queue_bef ore_assly_mc [assly_mc_no] . 

no_of_subasslys_into_assly ) ; 

f inished_assly__no: = 
aam[assly_mc__no] . 

f inished_assly_.no ; 

(* queue_before_assly_mc[assly_mc_no] . 

no_of_subasslys_into_assly: = 
abda[ finished_assly_.no] . 

no_o f _ch i 1 d_a s s 1 y s ; * ) 

for subassly_.no : =1 to max_subassly_no_Into_an_assly do 
begin 

queue_before_assly_mc [assly_mc_no] . 

details_of_queues [ subass ly_no] . queue_length: =0; 



quantity_available [subassly_no] : =0; (* added *) 

for quelg_no:=l to max_queue_length do 
begin 

queue_bef ore_assly_mc [assly_mc_no] . 

details_of_queues [subassly_no] . 

detail s_of _queue [ que 1 g_no ] . lot_no:=0; 
queue_bef ore_assly_mc [assly_mc_no] . 

detai ls_of_queues [subassly_no] . 

details_of_queue [quelg_no] . assly_no: =0 

end ; 
end; 
end; 

(* now all me queues and assly me queues are initialised *) 

(* now we do all pointer initialisations *) 

event_freetop: =nil; 
ptr_to_current_executable_event : =nil ; 
top_of_event_list: =nil; 
p: =nil ; 

(* pointer initialisations are complete *) 

(* initialising general machine data *) 

for mc_no:=l to max_mc_no do 
with agm[mc_no] do 
begin (* with *) 
tpt: =0. 0; 

total_idle_time: =0.0; 
cu: =0. 0; 
end; (* with *) 

(* end of initialising general machine data *) 

(* initialising assembly machine data *) 

for assly_mc_no: =16 to 20 do 
with aam[assly_mc_no] do 
begin (* with *) 
tpt : =0. 0; 

total_idle_time: =0. 0; 
cu: =0. 0; 
end; (* with *) 

(* end of initialising assembly machine data *) 

(* lot size needed or not is determined *) 

for assly_no:=l to max_assly_.no do 
abda[assly_no] . 

lot_sizing_to_be_done: =true; 
for assly_no:=l to max_assly_no do 
begin (* for *) 

if (abda[assly_no] . 

immediately_needs_processing) 



then 

begin (* then *) 

assly__no_af ter jprocessing: = 
abda[assly_no] . 

assly_no_af ter_processing; 
abda [assly__no_af terjprocessing] . 

lot_sizing_to_be__done: =false; 

end; (* then *) 
end; (* for *) 

for assly__no: =1 to max__assly_.no do 
alatasslyjno] . 

lot_no_which_will_enter_the_shop: =1 ; 

for mc_no : =1 to max__no_of jproces_on_any_assly do 

mc__busy [mc_no] :=false; 
end; (* initialisations *) 

procedure f ind__the_next__executable__event ; 

begin (* f ind_the__next_executable_event *) 

ptr_to_current_executable__event :=top_of_event_list; 
if (ptr_to_current_executable_event=nil) then 
flag: =true 
(* no event left *) 
else 
begin 

top__of_event__list : =top_of_event_list" . next__event; 
clock_time: =ptr_to_current__executable__event / ' . 

execution__time; 

end; 

end; (* f ind__the_next_executable_event *) 

procedure search__and_insert_f or__general_mc_queue (mc_no : 

var i , j , k, n, index, 

insert__position, 
assly__no , 
lot__no, 

cur rent_ass ly_.no, 
current__lot_no , 
next_assly__no, 
assly_mc__no : integer; 

priori ty_value_of_current_lot : real ; 

insert_position__found : boolean; 

begin (* search_and__insert_for__general_mc_,queue *) 

current_assly_no : =ptr_to_current_executable__event" . 

array_of _relevant_nos [ 1 ] ; 

current_lot__.no : =ptr__to_current_executable__event" . 

array_of_relevant__nos [2] ; 

if (queue__before__general_mc [mc_no] . queue__length = 0) 

then begin 

queue_bef ore__general_mc [mc__no] . 

details__of_queue [ 1 ] . assly_.no : = 
current_assly_.no ; 


integer) 
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queue_bef ore_general_mc [mc_no ] . 

detail s_of_queue [ 1 ] . lot_no: = 
current_lot_no ; 

queue_before_general_mc [mc_no] . queue_length: = 

queue_before_general_mc[mc_no] ,queue_length + 1; 

end 

else 

begin 

1 : =1 ; j : =queue_bef ore_general_mc [mc_no] . queue_length; 

current_assly_no: =ptr_to_current_executable_event~. 

array_of_relevant_nos [ 1 ] ; 

current_lot_no : =ptr_to_current_executable_event~ . 

array_of _re levant_nos [ 2 ] ; 

next_assly_no: =abda[current_assly_no] . assly_no_af ter_processing; 
assly_mc_no: =abda [next_assly_no] . assly_on_mc_no; 
if not first_run then 

coordinat ion_index [ cur rent_ass ly_no ] : = 

1 + ( latest_assly_t ime [ current_assly_no , assly_mc_no ] - 
prev_mean_time[ assly_mc_.no] )*0. 1; 
priority_value_of_current_lot: = 
ala [current_assly_no] . 

dol [current_lot_no] . priority_value * 

coordinat ion_index [ current_assly_no ] ; 
insert_position_found: =false; 

repeat 

k: =(i+j) div 2; 
assly_no: = 

queue_before_general_mc[mc_no] . details_of_queue[k] .assly_no; 
lot_no: = 

queue before g eneral me [me no] . details_of_queue[k] . lot_no; 

if (priority_value_of_current_lot > 
ala[assly_no] . 

dol [lot_no] .priority_value) 
then begin (* then *) 
if (k=l) 
then begin 

insert_position_found: =true; 
insert_position: =1 ; 
end 

else 

begin 

assly_no: = 

queue_before_general_mc [mc_no] . 

details_of_queue[k-l] . assly_no; 

lot_no: = 

queue before g eneral me [me no] . 

details_of_queue[k-l ] . lot_no; 
if (priority_value_of_current_lot <= 
ala [assly_.no] . 

dol [lot_.no ] ,priority_value) 

then begin 

insert_position_found: =true; 
insert_position: =k; 

end 

else j : =k-l ; 
end; 
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end (* then *) 
else begin (* else *) 

if (k-queue__bef or e__general__mc [mc__no] . queue length) 
then begin “ 

insert_position_found: =true; 
insert_position: = 

queue_before_general_mc[mc_no] . 

queue_length + 1 ; 

end 

else 

begin 

asslyjno : = 

queue_before_general_mc[mc_no] . 

detai ls_of__queue [k+1 ] . asslyjno ; 

lot_no:= 

queuejDefore_generaljnc[mc_no] . 

details_of_queue[k+l] . lotjno; 
if (priori ty_value_of_current_lot >= 
ala [asslyjno] . 

dol [ lot_.no] . priori ty__value) 

then begin 

insertjposition_found: =true; 
insert__position: =k+l ; 

end 

else i:=k+l; 
end ; 

end; (* else *) 
until ( insert jpos it ionjfound) ; 

for index: squeue_before_general_mctmc_.no] . queue_length 

downto 

insert j>osition 
do 

queue_bef ore_general jnc [mc_no ] . 

detai ls_of_queue [ index+1 ] : = 
queue_bef ore_general_mc [mc_no] . 

detail s_of _queue [ index] ; 
queuejbef ore_general_mc [mc_no ] . 

details_of_queue [ insert jsosition] . asslyjno: = 
current_assly_no ; 
queue_before__general jnc [mc_no] . 

detai ls_of_queue[insert_position] . lot jio: = 
current_lo t jno ; 

queue_before_general_mc [mcjno] . queue__length: = 
queue_before_general_mc[mc_no] . queue_length + 1; 
end; 

end; (* search_and_insert_f or_general_mc_queue *) 

procedure search_and_insert_f or_assemblyjnc_queue (mcjno : integer ) ; 

var f inish_assly_no , 
current__asslyjio , 
current_lot_no, 
que_len, 

u PP er _sear ch_l im , 
searehjio ; integer; 



positionjfound : boolean; 

begin (* search_and_insert_f or_assembly_jnc_queue *) 

current__assly_no : =ptr_to_current__executable_event /v . 

array_of_relevant_nos [ 1 ] ; 
current_lot_no: =ptr_to_current_executable_event~ . 

array__of_relevant_nos [ 2 ] ; 
f inish_assly_.no : =aam [mc_no ] . f inished_assly_no ; 

upper__search_lim: =abda[f inish__assly_no] . 

no__of_children_for_assembly; 
search_no: =1 ; posit ionjf ound: =false; 

while (search__no <= upper_search_lim) and not (position_found) do 
if (abda[f inish_assly_no] . 

arra y__°f- c hild_assly [search_no] . assly_no = 
current_assly_no ) 
then position_found: =true 
else search_no: =search_no + 1; 

if not (positionjfound) 
then begin 

wri teln( * asslyjio not found in insert queue for assly me’) 
halt; 
end 

else (* position has been found *) 
begin (* else *) 

queue_before_assly_mc [mc_no] . details__of ^queues [search_no] . 
queue_length: = 

queue_before_assly_mc [mc_no] . details__of_queues [search_no] . 

queue__length + 1 ; 
que_len: = 

queue_before_assly_mc [mc_no] . details_of_queues [search_no] . 
queue_length; 

queue_before_assly_mc [mejio] . details_of_queues [search_no] . 

details__of_queue [que_len] . assly_no: =current_assly_no; 
queue_before_assly_mc[mc_no] . details_of_queues [search_no] . 

details_of_queue [que_len] . lot_no: =current_lot_no; 
queue_before_assly_mc[mc_no] . 

quant ity__avail able [search_no] : = 
queue_before_assly_mc [mc_no] . 

quant ity__available [sear ch_no] + 
ala[current_assly_no] . 

dol [current_lot_no] . 
original_size; 
end; (* else *) 

end; (* search_and_insert_f or_assembly_mc_queue *) 

procedure execute_collect_statistics_from_now_on; 

var machine jno , 
queue jio , 
n °__pf ^queues , 
queue_length, 
queue_length_index , 
asslyjio, 



lot no 


: integer; 


(* collect_statistics_from_time *) 

begin (* execute_collect_statistics__f rom_now_on *) 

for machine_no: =1 to max_mc_no do 
begin (* for general__machine_no *) 
agm[machine__no] . tpt : =0.0; 
end; (* for general jmachine_no *) 
for machine_no: =16 to 20 do 
begin (* for machine_no *) 
aam [machine_no ] . tpt : =0 . 0 ; 
end; (* for machine_no *) 
for machine_no: =1 to max_mc_no do 
begin (* for general__machine_no *) 

queue_length: =queue_bef ore__general_mc [machine_no] . queue_length; 
for queue__length_index : =1 to queue_length do 
begin (* for queue__.no *) 

assly_no: =queue_before_general_mc[machine_no] . 

details__pf_queue [queue_length__index] . assly_.no 
lot_no : =queuejDefore__general__mc[machine_no] . 

details__of_queue [queue_length__index] . lot_no; 

ala[assly__no] . 

dol [lotjno] . total_wait_time_in_queues_so_far : =0. 0; 
ala[assly_no] . 

dol[lot__no] . time_of_entry__into_the_shop: = 
collect_statistics_from_time; 
end; (* for queue_.no *) 
end; (* for general_machine_no *) 
for machine__no: =16 to 20 do 
begin (* for assly_machinejno *) 

no__of__queues: =queue_bef ore_assly_mc [ machine_.no ] . 

no__of__subasslys_into_assly ; 
for queue_no:=l to no_of_queues do 
begin (* for queue__no *) 

queue_length: =queue_bef ore_assly_mc [machine_.no] . 

details_of_queues [queue__no] . queue_length; 
for queue_length_index: =1 to queue_length do 
begin (* for queue_length__index *) 
assly — no: = 

queue_bef ore_assly_mc [machine_no ] . 
details_of_queues [queue_.no ] . 

details_of_queue [queue_length_index] . assly_no; 
lotjno: = 

queue_bef ore_assly_mc [machine_no] . 
details_of_queues [queue_no] . 

detai ls_of_queue [queue_length_index] . lot_no; 
ala [assly_no] . 

do 1 [ lot no ] . to tal__wai t_t ime__in__queues_so_f ar : =0 , 0 ; 

ala [assly_.no] . 

dol [ lot_no] . time__of__entry__into_the._shop: = 
collect_statistics_jfrom_time; 
end; (* for queue_length__index *) 
end; (* for queue_.no *) 
end; (* for assly_machine_.no *) 
end; (* execute__collect__statistics_jfrom._now_pn *) 
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procedure execute_lo t_enters_the_shop ; 

var assly_no, 
mc_no , 

lot_no : integer; 

begin (* execute_lot_enters_the_shop *) 

newmodif ied (new_event ) ; 

new_event~ . array_of_relevant_nos [ 1 ] : = 

ptr_to_current_executable_event" . array_of_relevant_nos [ 1 ] ; 
new___event~ . array_of_relevant_nos [2] : = 

ptr_to_current_executable__event^ . array_of_relevant_nos [2] ; 
new_event~ . execution_time: = 

ptr_to__current_executable__event'\ execution_time; 

assly_no : =ptr_to_current_executable_event >N . array_of_relevant_nos [ 1 ] ; 
1 o t_no : =p t r__to_current_executab 1 e_event * . ar ray_of _r e 1 evant_nos [ 2 ]; 
if (abda[assly__no] . 

immediate ly__gets_into__assly) 

then begin 

new_event~ . event : =lot_enters__the_assly__mc__queue ; 
mc_no: =abda [assly_no] . 

assly_on_mc_no ; 

new_event" . array_of_relevant_nos [3] : =mc__no; 
end 

else if (abda [assly_no] . 

immediately_needs_processing) 
then begin 

new__event". event : =lot_enters_the__general_mc__queue ; 
mc__no: =abda [asslyjno] . 

aps[l] . mcjio; 

new_event^ . array__of__relevant jnos [3] : =mc_no; 
end 

else new_event~ . event : =lot_enters_the_f inished_product_list ; 
ala[assly_no] . 

dol [lot_no] . time_of_entry_into_the_shop 

: =ptr__to_current_executable_event" . execution_time; 
plug_in_the_event (new_event ) ; 

end; (* execute_lot_enters_the_shop *) 

procedure execute_lot_enters_the__general_mc_queue; 

var assly_no, 
lot_no, 
queue__lg, 
mc_no__to_visit , 
process_no : integer; 

begin (* execute_lot_enters_the_general_mc_queue * ) 

assly_no : = ptr_to_current_executable_event~ . array__of_relevant__nos [13; 
lot___no: =ptr_to_current___executable_event ^ . array__of__relevant_nos [2] ; 
mc_jno___to_visit : =ptr_ to___current_executable__event A . 

array_of_relevant_nos[3] ; 

process__no: =ala [assly_no] . dol [ lot_no] . no_of_agm_visited + 1; 
ala[assly_no] . 

dol [lot__no] . current_mc__no: =mc_no_to_visit ; 



ala[assly_no] . 

dol [lot_no] . time_of_joining_the_queue: = 
ptr_to_current_executable_event* . execution_time; 

ala[assly_no] . 

dol [lot_no] . cpt : = 
ala [assly_.no ] . 

dol [lot_.no] . original_size * 
abda [assly_.no ] . 
aps [processjno] . 

process ing_time_per_unit 
+ 

abda [ass ly_no] . 
aps [process_.no] . 
setup_time; 

search_and_insert_for_general_mc_queue(mc_no_to_visit) ; 
queue_lg: =queue_before_general_mc[mc_no_to_visit] . queue_length; 
if (queue_lg = 1) and (mc_busy [mc_no_to_visit]=false) then 
begin 

newmodif ied (new_event ) ; 

new_event" . event: =processing_on_general_mc_starts; 
new_event" . array_of_relevant_nos [ 1 ] : =assly_no ; 
new_event" . array_of_relevant_nos [2] : =lot_no; 
new_event" . array_of_relevant_nos [3] : =mc_no_to_visit ; 
new_event~. execution^ ime: = 

ptr_to_current_executable_event" . execution_time; 
pl u g_in__the_event (new_event ) ; 
end; 

end; (* execute_lot_enters_the_general_mc_queue *) 

function process ing_can_s tar t_on_assly_mc 

(mc_no_to_visit : integer) : boolean; 

var processing_can_start : boolean; 

pos_of_lot_no , 

lot_size, 

no_of __chi ldren , 

chi ldren_index , 

quant i ty_to_be_cont r ibuted , 

f inish_assly_no : integer; 

begin (* process ing_can_s tar t_on_assly_mc *) 

finish_assly_.no := aam[mc_no_to_visit] . 

f inished_assly_no 

{queue_before_assly_mc [mc_no_to_visit ] . f inish_assly_no; > 
pos_of_lot_no : = 
ala[f inish_assly_no] . 

lo t_no_which_wi 1 l_enter_the_shop ; 
lot_size: = 

ala[f inish_assly_no] . 

dol [pos_of_lot_no] . original_size; 
no_of_chi ldren: =abda[f inish_assly_no] . 

no_o f _ch i 1 d r en_f o r_a s s emb 1 y ; 
processing_can_start : =true; 
children_index: =1 ; 
while (processing_can_start ) 
and 



(children_index <= no_of_children) 
do 

begin (* while *) 

quantity_to_be_contributed: = 
abda[f inish_assly_no] . 

array_of_chi ld_assly [chi ldren_index] . qty_to_be_contributed 
* lot_size; 

if (quantity_to_be_contributed > 

queue_before_assly_mc [mc_no_to_visit] . 
quant i ty_ava i 1 ab 1 e [ ch i ldren_index ] ) 
then processing_can_start: =false 
else children_index: =children_index + 1 ; 
end; (* while *) 

processing_can_start_on_assly_mc: =processing_can_start; 

end; (* processing_can_start_on_assly_mc *) 

procedure execute_lot_enters_the_assly_mc_queue ; 

var assly_no, 

f inished_assly_no , 
pos_of_lot_no, 
lot_no, 
queue_lg, 

mc_no_to_visit : integer; 

begin (* execute_lot_enters_the_assly_mc_queue *) 

assly_no : =ptr_to_current_executable_event - ' . array_of_relevant_nos [ 1 ] 
lot_.no : =ptr_to_current_executable_event~ . array_of_relevant_nos [2] ; 
mc_no_t o_v i s i t : =p t r_t o_cur r en t_exe cut ab 1 e_event ~ . 

array_of_relevant_nos[3] ; 
assly_count [mc_no_to_visit] : = 

assly_count [mc_no_to_visit] + 1; 
ala[assly_no] . 

dol[lot_no] . time_of_joining_the_queue: = 

ptr_to_current_executable_event A . execution_time; 
time_of_reaching[assly_no, lot_.no ] : = 

ptr_to_current_executable_event' s . execution_time; 
present_mean_time [mc_no_to_visit ] : = 

( (assly_count [mc_no_to_visit ] -1 ) *present_mean_time [mc_no_to_visit ] 
time_of_reaching [assly_no, lot_no] )/assly_count [mc_no_to_visit] ; 
if (time_of_reaching[assly_no, lot_no]> 

present_mean_time [mc_no_to_visit] ) 

then 

latest_assly_time[assl y_no , mc_no_t o_v i s i t ] : = 
time_of_reaching[assly_no, lot_no] ; 
arrival_time[times+l, assly_.no, lot_no] : = 

time_of_reaching[assly_no, lot_.no] ; 

sear ch_and_inse r t_f or_as semb 1 y_mc_queue ( mc_no_t o_v i s i t ) ; 

if (processing_can_start_on_assly_mc(mc_no_to_visit] ) 

then 

begin (* then fire the event processing_starts_on_assly_mc *) 
if (ptr_to_current_executable_event~ . execution_time< 

earliest_lot_time [mc_no_to_visit] ) then 
earliest_lot_time [mc_no_to_visit] : = 
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ptr_to_current__executable_event" . execution time 

else if ~ 

(ptr_to__current_executable__event^ . execution^ ime> 

latest__lot_time [mc_no_to_visit] ) then 
latest_lot__time [mcjio_to_visit ] : = 

ptrJ:o_current_executable_event~.execution_time; 

f inished_assly_jio: = 
abda[assly_no] .parent_no; 
newmodif led (new_e vent ) ; 

new_event" . event : =processing_on_assly_mc_starts; 
new_event A . array_of_relevant_nos [1 ] : =f fnished_assly_no; 
pos_of_lot__no : = 
ala [f inished_assly_no] . 

lotjno_which_will_enter_the_shop; 
new_event A . array_of__relevant_nos [2] : =pos_of_lotjno; 
new__event A . array_ofj~elevant_nos [3] : =mc_no_to_visit ; 
new_event" . execut longtime: = 

ptr_to_current_executable_event" . execution_time; 
plug_in_the_event (new_event ) ; 

end; (* then fire the event processing_starts_on_assly_mc *) 

end; (* execute_lot_enters_the__assly_mc_queue *) 

procedure execute_processing_on_general_mc_s tarts ; 

var assly__no, 
lotjno, 
queue_lg, 
queue_lg__no , 
size, 

mc_no__to_yisit : integer; 

processor inish_time : real; 

begin (* execute_j>rocessing_on_general__mc_starts *) 

assly_no : =ptr_to < _current_executable_event'\ array_of_relevant_nos [ 1 ] ; 
lot_no : =ptr_to_current_executable_event" . array_of_relevant_nos [2] ; 
mc_no__to__visi t : =ptr_to_current_executable_event" . 

array___of_relevant_nos[3] ; 
mc_busy [mc_no_to_visit ] : = true; 
size: =ala[assly__no] . 

dol [ lot__no] . 

original_size; 
processor inish_time: = 

ptr_to_current_executable_event A . execution__time + 
ala [assly_no] . 
dol [lot_.no] . cpt ; 

(* now insert an event called processing gets over on the same machine *) 
newmodif ied(new_event ) ; 

new_event~ . event: =processing_on_general i jnc_gets_over ; 
new -_event~ . array_ofjrelevant_nos [ 1 ] : =assly_no ; 
new_event~ . array_of_relevant_nos [2] : =lot_no; 
new__event" . array_of_relevant_nos [3] : =mc_no_to_visit ; 
new__event^ . execution^ ime: =process_f inish_time; 
plug_in_the._event (new_event ) ; 



(* the event has been inserted *) 


queue_lg:=queue_before_general_mc[mc_no_to_visit] . queue_length; 

for queue_lg_jio: =2 to queue_lg do 
queue_before_generaljnc [mcjno_to_visit ] . 

detailsjD fjqueue [queue_Jgjio-l ] : = 
queue Joefore^generaljmc [mcjno__toj/isit ] . 

details_of__queue [queue_lg_no] ; 
queue_before_generaljnc [mcjno_toj/isit ] . queue__length: = 
queueJbefore_^generaljrnc[mcjio_to_visit] . queue_length - 1; 

(* now we modify the assly-no and lot-no data *) 

alatasslyjno] . dol [lotjno] . time_whenjprocessingj3tarts: = 
ptr_to__current_executable_event~ . execut longtime; 
with alatasslyjno] . 

dol [lotjno] do 
begin (* with *) 

current jncjno: =mcjno_to_visit ; 
total j^ait_time_injqueues_soJf ar : = 
total_wait_time_injqueuesj3ojfar + 

{ ptr_to_current_executable_event* . executionjime -> 
timejwhenjprocess ingest arts - 

t ime_of _ j o ining the jqueue ; 

end; (* with * ) 

end; (* executejprocessing_on__general jnc_starts *) 

procedure executejprocessing_on_asslyjnc_starts; 

var asslyjno, 
lotjno, 
queue_lg, 
size, 
lot_size, 

last_quantity__available, 

quant i t y_ava i 1 ab 1 e , 

quant i ty__to_be_cont r ibuted , 

children_index, 

n°__of __chi ldren , 

lowlim, 

lowlim_index, 

1 o t_no__o f _ch i 1 d , 
asslyjno__of _chi Id , 
q ty jp i cked_up , 

mc_no_toj/isit : integer; 

processjf inishjtime : real; 

begin (* execute jprocessingj>n_asslyjnc_starts *) 

asslyjno : =ptr_to_current_executable_event^ . array_of__relevant__nos [ 1 ] 
lotjno : =ptr_to__current__executable__event~ . array_of jrelevant jios [2] ; 
mcjno_to_visi t : =ptr_to_current_executable_event" . 

arrayjnf__relevant jios [3] ; 

lot_size: = 
ala [asslyjno] . 

dol [lotjno] . original jsize; 
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no_of_children: =abda[assly_no] . 

no _° f _ch i 1 d r en_f o r_as semb 1 y ; 
for children_index:=l to no_of_chi ldren do 
begin (* for children_index *) 
quant i ty_to_be_contributed : = 
abda[assly_no] . 

array_of_child_assly [children_index] . qty_to_be_contributed 
* lot_size; 
lowlim: =1 ; 
lo t__no_of_chi Id : = 

queue_before__assly_jnc [mc_no_to_yisit ] . 
detail s_pf_queues [children_index] . 
details___of_queue [lowlim] . 
lotjno; 
asslyjno_of__child: = 

queue__bef ore_assly__mc [mc_no__to_yisit ] . 
details_of_queues [children_index] . 
details_of_queue [ lowlim] . 
assly__no; 

last__quant i ty__avai lable : =0 ; 
quanti ty__avai lable : = 
ala[assly_no_of_child] . 

dol [ lot jno_of__child] . qty_left; 
while (quantity_available < quant ity_tojDe_contributed) do 
begin (* while *) 

qtyjpicked_up: ==quantity_available-last_quantity_avai lable; 
ala[assly_no_of_child] . 

dol [ lot_no_of_child] . 

to tal_wai t_t ime_in_queues_so_f ar : = 
ala [ assly__no__of__chi Id ] . 

dol [ lot_no_of_child] . 

total_wai t_time_in_ > queues_sojfar + 

( (ptr_to_current_executable_event~ . execution^ ime - 
ala[asslyjrio_pf_child] . 

dol [ lotjno_of_child] . 

t ime_of_entry_into_the_shop)*qty_picked_up/ 
ala[assly_no__of_child] . 

dol [lot_no_of_child] . original_size) ; 
lowlim: =lowlim + 1; 
lot_no__of__child: = 

queue_before_assly_mc [mc_no_to_visit ] . 
detail s_of_queues [childrenJLndex] . 
cietai ls__of__queue [ lowlim] . 
lot_no; 

asslyjno_of__child: = 

queue Jbef or e_assly__mc [mc_no_to_visit] . 
details_of_queues [childrenJLndex] . 
details_of_queue [ lowlim] . 
asslyjno; 

1 as t_quan t i t y_a va i lable: =quant i t y_ava i lable; 
quant ity_avai lable : = 
quant ity_ava liable + 
ala[asslyjno_of_child] . 

dol [lot_no__of_child] . qty_left; 
end; (* while *) 

if ( quant ity_ava liable = quant ity_to_be_contributed) 

then 

else (* quant ity_ava liable > quantity_to_be_contributed *) 



begin 

qty_picked_up:=quantity_available-last_quantity_available 

ala[assly_no_of_child] . 

dol [lot_no_of__child] . qty_left:= 
quant i ty__avai 1 able - quant ity_to___be__contributed; 
ala[assly_no_of__child] . 
dol [ lot_no_of_child] . 

to tal_wai t__t iine__in_queues_so_f ar : = 
ala[assly__no__of__child] . 
dol [ lot_no_of__child] . 

total_wait__t ime__in_queues__so_JTar + 

( (ptr_to_current_executable_event~ . execution__time - 
ala[assly_no_of_child] . 
dol [lot_no_of_child] . 

time_of_entry_into_the_shop)*qtyjpicked_up/ 
ala[assly_no_jDf_child] . 

dol [ lot_no_of_child] . original__size) ; 
lowlim: =lowlim - 1; 
end; 

(* lowlim represents now the nos. in queue which get absorbed 
in the assly *) 

queue_lg: = 

q u e u e_b e f o r e__a s s .1 y __m c [mc__no__to__visit 3 . 
detai ls_ofjqueues [children_index] . 
queue_length; 

for lowlim_index: = (lowlim+l ) to queue_lg do 
queuejbef ore_asslyjmc [mc_no_to_visi t ] . 
detai ls_ofjqueues [chi ldren_index] . 

detai ls__of_queue [ lowlim_index-lowlim] : = 
queue__bef ore_assly_mc [mc_no_to___visit ] . 
details___of_queues [children__index] . 
detai ls_of_queue [ lowlim_index] ; 
if (queue__lg >= lowlim) then 
queuejbef ore_asslyjnc [mcjno_to_visit ] . 
detai ls__of_queues [chi ldren_index] . 
queue__length: =queue__lg - lowlim 

else 

queuejbef ore_asslyjmc [mc_no_to_visit ] . 
detai ls_of_jqueues [chi ldren_index] . 
queue_length: =0; 

queue__bef ore_assly_mc [mc_no__to_visi t ] . 
quant ity_available [children_index] : = 
quant ity_avail able - quant ity_to_be_contributed; 
end; (* for childrenJLndex *) 

processJF inish_time: = 

ptr_to_current_executable_event / " . executionjime + 
lot_size * 

aam [ mc_no__ t o_v isit] .process! ng_t i me_p e r_un i t ; 
newmodif ied (new_event ) ; 

new_event" . event : =processing_on_assly_mc_gets_over; 
new^event^ . array_ofjrelevantjnos [ 1 ] : =assly_no; 
new^event^ . array_ofjrelevant_nos [ 2 ] : =lot_no; 
new^event"'' . array_of_relevant_nos [3] : =mc_no_to__visit ; 
new_event" . execution_time: =process_f inish_time; 



plug__in_the_event (new__event ) ; 

end; (* execute_processing__on_assly_mc_starts *) 
procedure execute_jprocessing_.on_general_mc_gets_pver ; 

var assly_.no , 
lot__no, 
queue_lg, 
queue_lg_.no, 
size, 

next_process_no , 
mc_no_to_visit : integer; 

begin (* execute_processing_on_general_mc_gets_over *) 

assly_.no : =ptr_to_current_executable_event~ . array_of_relevant_nos [ 1 ] 
lot_.no : =ptr_to_current_executable_event~ . array_of_relevant_nos [2] ; 
mc_no_to_visit : =ptr_to_current_executable_event" . 

array_of_relevant_nos [3] ; 
mc_ busy [mc_no_to_visit ] :=false; 

if (queue_before_general_mc [mc_no_to_visit ] . queue_length > 0) then 
begin 

newmodif ied(new_event ) ; 

new_event " . event : =processing_on_general_mc_starts ; 
new_event" . array_of_relevant_nos [1] : = 

queue_bef o r e_gene r a l_mc [ mc_no_to_v i s i t ] . 

details_of_queue[l] . assly_.no ; 
new_event" . array_of_relevant_nos [2] : = 

queue_before_general_mc [mc_no_to_visit ] . 

details_of_queue[l] . lot_.no; 
new_event" . array_of_relevant_nos [3] : =mc_no_to_visit ; 
new_event" . execution_time: = 

ptr_to_current_executable_event~ . execution_time; 
plug_in__the_event (new__event ) ; 
end; 

(* beginning of initialisations and updations *) 
ala[assly_no] . 

dol[lot_no] . no_of_agm_visited 
ala [assly_.no] . 

dol [lot_no] . no_of_agm_visited + 1; 
ala [assly_.no] . 

dol [lot_.no] . tpt_completed_so_far : = 
ala[assly_no] . 

dol[lot_no] . tpt_completed_so_far + 
ala [assly_.no ] . 

dol [lot_.no] . cpt ; 
agm[mc_no_to_visit ] . tpt : = 
agm[mc_no_to_visit] .tpt + 

ala [assly_no] . 
dol [lot_.no] . cpt; 


ala [assly_no] . 

dol [lot_.no] . priority_value: = 
(ala[assly_no] . 



dol [lot_no] . total_required_processing_time - 

ala [assly_no] . 

dol[lot_no] . tpt_completed_so_far - 

(clock_time - 
ala [assly_no] . 

dolt lot_.no] . required_completion_time) )/ 
ala[assly_no] . 

dol [lot_no] . total_required_processing_time; 

(* end of Initialisations and updations *) 

if (ala[assly_no] . 

dol [ lot_no] . 

no_o f agm_v i s i t e d ) = 

(abda [assly_no] . 

no_of jprocessesjneeded ) 
then (* all the processes are done *) 
begin (* then *) 
ala [assly_.no ] . 
dol [ lot_.no ] . 

process ing_complete: =true; 
ala [assly_no] . 

dol [ lot_no ] . actual_completion_time: = 
ptr_to_current_executable_event" . execution_time; 
assly_no: = 
abda [ ass ly_no] . 

assly_no_af ter_processing; 
ala[assly_no] . 

dol[lot_no] . time_of_entry_into_the_shop: = 
ptr_to_current_executable_event~ . execution_time; 
if (abda [ass ly_no] . parent_no = 0) 
then (* joins the finished goods queue *) 
begin 

newmodif ied (new_event ) ; 

new_event" . event : =lot_enters_the_f inished_product_list ; 
new_event" . array_of_relevant_nos [ 1 ] : =assly_no; 
new_event~ . array_of_relevant_nos [2] : =lot_no; 
new_event" . array_of_relevant_nos [3] : =0; 
new_event" . execution_time: = 

ptr_to_current_executable_event" . execution_time; 
plug_in_the_event (new_event ) ; 
end 

else (* joins the assly machine queue *) 
begin 

mc_no_to_visit : = 

abda [ass ly_no] . assly_on_mc_no; 

newmodif ied (new_event ) ; 

new_event" . event: =lot_enters_the_assly_mc_queue; 
new_event " . array_of_relevant_nos [ 1 ] : =assly_no ; 
new_event~ . array_of_relevant_nos [2] : =lot_no; 
new_event" . array_of_relevant_nos [3] : =mc_no_to_visit ; 
new_event" . execution_time : = 

ptr_to_current_executable_event^ . execution_time; 
plug_in_the_event (new_event ) ; 
end; 

end (* then *) 

else (* joins the next machine *) 
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if (ala[assly_no] . 

dol [ lot_no] . 

no_of_agm_visited) < 

(abda [assly_no] . 

no_of_processes_needed ) 

then 

begin (* else *) 
next_process_no : = 
ala[assly_no] . 

dol [ lot_no] . 

no_of_agm_visited+l ; 
mc_no_t o_v i s i t : = 
abda [ ass ly_no] . 

aps [next_process_no] .mc_no; 
newmodif ied (new_event ) ; 

new_event'‘ . event : =lot_enters_the_general_mc_queue; 
new_event" . array_of_relevant_nos [ 1 ] : =assly_no; 
new_event~ . array_of_relevant_nos [2] : =lot_no; 
new_event~ . array_of_relevant_nos [3] : =mc_no_to_visit ; 
new_event"' . execution_time : = 

ptr_to_current_executable_event~ . execution_time; 
plug_in_the_event (new_event ) ; 
end; (* else *) 

end; (* execute processing on general me gets_over *) 

procedure execute_processing_on_assly_mc_gets_over ; 

var asslyjno, 

f inished_assly_no , 
lot_no, 

next_jprocess_no , 
mc_no_to_visit : integer; 

begin (* execute_processing_on_assly_mc_gets_over *) 

assly_no : =ptr_to_current_executable_event - ' . array_of_relevant_nos [ 1 ] ; 
lot_no: =ptr_to_current_executable_event"' . array_of_relevant_nos [2] ; 
mc_no_to_visit : =ptr_to_current_executable_event~ . 

array_of_relevant_nos [3] ; 

ala[assly_no] . 

lo t_no_whi ch_wi 1 l_enter_the_shop : = 
ala[assly_no] . 

lot_no_which_will_enter_the_shop + 1; 
ala [assly_no ] . 

dol[lot_no] . time_of_entry_into_the_shop: = 
ptr_to_current_executable_event / ' . execution_time; 
ala[assly_no] .dol [lot_no] . actual_completion_time: = 
ptr_to_current_executable_event / ' . execution_time; 

if ( abda [ ass ly_no ] . 

immediately g ets into.assly) 

then 

begin (* gets into assly *) 
f inished_assly_no: = 
abda [ ass ly_no] .parent_no; 
mc_no_to_v i s i t : = 

abda [f inished_assly_no] . assly_on_mc_no; 
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newmodif ied (new_event ) ; 

new_event~ . event : =lot_enters_the_assly_mc_queue; 
new_event~ . array_of_relevant_nos [ 1 ] :=assly_no; 
new_event~ . array_of_relevant_nos [2] : =lot_no; 
new_event~ . array_of_relevant_nos[3] : =mc_no_to_visit; 
new_event~ . execution_time: = 
ptr_to_current_executable_event'' . execution_time; 
plug_in_the_event (new_event ) ; 
end (* gets into assly *) 
else if (abda [assly_no] . 

i immediate ly_needs_processing ) 

then 

begin (* processing of the assembly starts *) 
mc_no_to_visit : = 
abda [ass ly_no] . 

aps [ 1 ] . mc_no ; 
newmodif ied (new_event ) ; 

new_event"' . event : =lot_enters_the_general_mc_queue; 
new_event ,N . array_of_relevant_nos [ 1 ] : =assly_no; 
new_event~ . array_of_relevant_nos [2] : =lot_no ; 
new_event~ . array_of_relevant_nos [3] : =mc_no_to_visit ; 
new_event~. execution_time:= 

ptr_to_current_executable_event'\ execution_time; 
plug_in_the_event (new_event ) ; 
end (* processing of the assembly starts *) 
else 

begin (* gets into the finished product list *) 
newmodif ied (new_event ) ; 

new_event~. event :=lot_enters_the_finished_product_l 1st; 
new_event" . array_of_relevant_nos[l] : =assly_no; 
new_event~ . array_of_relevant_nos [2] : =lot_no; 
new_event'~ . array_of_relevant_nos [3] : =0; 
new_evenf' . execution_time: = 

ptr_to_current_executable_event" . execution_time; 
plug_in_the_event (new_event ) ; 
end; (* gets into the finished product list *) 

end; (* execute_processing_on_assly_mc_gets_over *) 

procedure execute_the_event ; 

begin (* execute_the_event *) 

case ptr_to_current_executable_event". event of 

collect_statistics_f rom_now_on : 

execute_collect_statistics_from_now_on; 
lot_enters_the_shop : 

execute_lot_enters_the_shop; 
lot enters the g eneral me queue : 

execute_lot_enters_the_general_mc_queue; 
lot_enters_the_assly_mc_queue : 

execute_lot_enters_the_assly_mc_queue; 
processing_on_general_mc_s tarts : 

execute_processing_on_general_mc_starts; 

processing_on_assly_mc_starts : 

executejprocessing_on_assly_mc_starts; 

processing on general me gets_over : 
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execute_processing_on_general_mc_gets_over ; 

processing_on_assly_mc_gets_over : 

execute_processing_on_assly_mc_gets_over; 
lot_enters_the_f inished_product_list . 

execute_lot_enters_the_finished_product list- 

end; (* case *) 

event_addf ree 1 i s t (p t r_to_cur r ent_executab le_event ) ; 

end; (* execute_the_event *) 

procedure co llect_and_print_s tat i sties; 

var assly_no, 
lot_no, 

no_of_orders_placed, 
order_no , 

current_demand_qty, 

current_demand_per iod , 

current_supply_qty, 

current_supply_period, 

machine_no, 

queue_length, 

queue_length_index, 

no_of _queues , 

queue_no , 

no_of_lots , 

general_mc , 

qty, period , 

mc_.no : integer; 

setup_cost , 

total_setup_cost , 

setup_cost_reductions, 

total_aggregate_lateness , 

inventory_carrying_cost , 

total theoretical inventory carryin g cost, 
theoretical_inventory_carrying_cost , 
total_inventory_carrying_cost : real; 
sum_util , 
av_ut i 1 , 
sum_range , 

coordination_range :real; 

range: array [range_of_assly_mcs] of real; 

begin (* collect_and_print_statistics *) 

(* calculating theoretical inventory carrying costs *) 

theoretical_inventory_carrying_cost : =0. 0; 
for assly_no:=l to max_assly_no do 
begin (* for assly_.no *) 

(* theoretical_inventory_carrying_cost: = 
theoretical_inventory_carrying_cost + 
array_for_output_of_w_w[assly_no] . holding_cost; 
theoretical_cost [assly_.no] : = 

arra y_f or_output_of_w_w[assly_no] . holding_cost; * ) 

qty:=o ; 
period: =1 ; 

for order_no:= 1 to array_for_output_of_w_w[assly_no] . 
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no _of_orders_placed do 

begin 

qty : =qty+array_for_output_of_w_w [ ass ly_no ] . 

order_schedule[order_no] . qty_produced- 
while (qty-demand[assly_no, period] >0) and 

(period<max_period_no) do 

begin 

theoretical_cost [assly_no] : = 

(qty-demand [assly_no, period] ) * 

abda[assly_no] . cost*inventory_cost_fraction 

+ theoretical_cost [assly_no] ; 
theoretical_inventory_carrying_cost: = 
theoretical_inventory_carrying_cost + 

(qty-demand [assly_no, period] ) * 

abda [assly_no] . cost*inventory_cost_fraction; 

period: =period+l ; 

qty: =qty-demand[assly_no, period] ; 
end; (* while qty >*) 
end; (* for- order_no *) 
no_of_lots: =ala [assly_no] . no_of_lots; 
for lot_no:=l to no_of_lots do 
begin (* for lot_no *) 

actual_cost [assly_no] : =actual_cost (assly_no] + 
ala[assly_no] . 

dol [ lot_.no ] . tot.al_wait_time_in_queues_so_far * 
abda [assly_no] . cost * inventory_cost_f raction ; 
end; (* for lot no * ) 
end; (* for assly_no *) 

(* end of calculating theoretical inventory carrying costs *) 

(* adding waiting times for asslys in queues also *) 

for machine_no: =1 to max_mc_no do 
begin (* for general_machine_.no *) 

queue_length: =queue_before_general_mc[machine_no] . queue_length; 
for queue_length_index: =1 to queue_length do 
begin (* for queue_no *) 

assly no:=queue before g eneral me [machine_no] . 

details_of_queue [queue_length_index] . assly_no; 
lot_no: =queue_before_general_mc[machine_no] . 

details_of_queue [queue_length_index] . lot_no; 

ala[assly_no] . 

dol [ lot_no] . total_wait_time_in_queues_so_far: = 
ala[assly_no] . 

dol [lot_no] . total_wait_time_in_queues_so_far + 
simulation_time - 
ala[assly_no] . 

dol[lot_no] . time_of_joining_the_queue; 
end; (* for queue_no * ) 
end; (* for general_machine_no *) 
for machine_no: =16 to 20 do 
begin (* for assly_machine_no *) 

no_of_queues: =queue_before_assly_mc [machine_no] . 

no_of_subasslys_into_assly; 

for queue_no:=l to no_of_queues do 
begin (* for queue_no *) 

queue_length: =queue_before_assly_mc[machine_no] . 



details_of_queues [queue_no] . queue_length - 
for queue_length_index: =1 to queue_length do 
begin (* for queue_length_index *) 
assly_no: = 

queue_before_assly_mc [machine_no] . 
details_of_queues [queue_no] . 

details_of_queue[queue_length_index].assly_no; 

lot_no:= 

queue_before_assly_mc [machine_no] . 
details_of_queues [queue_no] . 

details_of_queue [queue_length_index] . lot_no; 

ala [assly_no] . 

dol [lot_no] . total_wait_time_in_queues so far: = 
ala[assly_no] . 

dol [lot_no] . total_wait_time_in_queues_so_far + 
simulation_time - 
ala [assly_no] . 

dol[lot_no] . time_of_joining_the_queue; 
end; (* for queue_length_index *) 
end; (* for queuejno *) 
end; (* for assly_machine_no *) 

(* end of adding waiting times for asslys in queues also *) 

(* calculating actual inventory carrying costs *) 

total_inventory_carrying_cost : =0.0; 
for assly_.no :=1 to max_assly_no do 
begin (* for assly_.no *) 

no_of_lots : =ala[assly_no] . no_of_lots; 
for lot_no:=l to no_of_lots do 
begin (* for lot_no * ) 

inventory_carrying_cost: = 
ala[assly_no] . 

dol [lot_no] . total_wait_time_in_queues_so_far * 
abda [assly_no] . cost * inventory_cost_fraction/ 

no_of_hours_in_a_period; 

{ inventory_carrying_cost/no_of_hours_in_a_year; } 
ala[assly_no] . 

dol [ lot_no] . actual_inventory_carrying_cost: = 
inventory_carrying_cost ; 
total_inventory_carrying_cost : = 

total_inventory_carrying_cost + inventory_carrying_cost; 
abda [assly_.no ] . 

actual_inventory_carrying_cost : = 
abda [assly_.no] . 

actual_inventory_carrying_cost + 
inventory_carrying_cost ; 
end; (* for lot_no *) 
end; (* for assly_no *) 

(* end of calculating actual inventory carrying costs *) 

(* calculating setup costs *) 
total_setup_cost : =0 . 0 ; 
setup_cost_reductions: =0.0; 
for assly_no:=l to max_assly_no do 
begin (* for *) 



setup_cost : =array_for_output_of_w_w[assly_no] . order_cost; 
total_setup_cost: =total_setup_cost + setup_cost; 

no_of_orders_placed: = ~ 

array_f or_output_of_w_w[assly_no] • no_of_orders_placed; 
for order_no:=l to no_of_orders_placed do 
begin (* for order_no *) 

if ( ar ray_f or_output_of _w_w [ass ly_no ] . 

order_schedule [order_no ] .period_no * 
no_of_hours_in_a_period < 
collect_statistics_from_time) 
then setup_cost_reductions: = 
setup_cost_reductions + 
abda[assly_no] . total_setup_cost 
else if (array_for_output_of_w_w[assly_no] . 

order_schedule [order_no] . period_no * 
no_of_hours_in_a_period > 
simulation_time) 
then setup_cost_reductions: = 
setup_cost_reductions + 
abda [assly_no] . 

total_setup_cost ; 
end; (* for order_no *) 
end; (* for *) 

total_setup_cost : =tot.al_setup_cost - setup_cost_reductions; 

(* end of calculating setup costs * ) 

(* calculating aggregate_lateness *) 

total_aggregate_lateness: =0. 0; 
for assly__no:=l to max_assly_no do 
begin (* for assly_no *) 

no_of_lots: =ala [assly_no] . no_of_lots; 
for lot_no:=l to no_of_lots do 
begin (* for lot_no *) 
if (assly_no = 1) then 

ala[l ] . dol [ lot_no] . actual_completion_time: = clock_time; 
ala [assly_no] . 

dol [ lot_no] . aggregate_lateness_of_lot : = 
(ala[assly_no] . 

dol [ lot_.no ] . actual_completion_time - 
ala[assly_no] . 

dol [lot_no] . required_completion__time) * 
ala [assly_no] . 

dol [lot_no] . original_size; 
total_aggregate_lateness : = 
total_aggregate_lateness + 
ala t assly_.no ] . 

dol [ lot_.no ] . aggregate_lateness_of_lot; 
end; (* for lot_no *) 
end; (* for assly_no *) 

(* end of calculating aggregate lateness *) 


(* calculating units_lateness *) 
total_uni ts_lateness : =0 . 0 ; 



for assly_.no: =1 to max_assly_no do 
begin (* for assly_no *) 

no_of_lots: =ala [assly_no ] . no_of_lots; 
current_demand_qty : =0 ; 
current_demand_period: =0; 
for lot_.no :=1 to no_of_lots do 
begin (* for lot_no *) 
current_supply_qty: = 
ala [assly_.no ] . 

dol[ lot_.no] . original_size; 
current_supply_period: = 
trunc(ala[assly_no] . 

dol [lot_.no ] . actual_completion_time/ 
no_of_hours_in_a__period) + 1; 
while (current_supply_qty > 0) 

and 

(current_demand_period < max_period_no) 
do 

begin (* while current_supply_qty > 0 *) 

current_demand_period: =current_demand_period + 1; 
cur ren t_demand_q ty : =demand [assly_no, current_demand_jperiod] ; 
if (current_supply_period > current_demand_period) 
then 

begin C* then *) 
ala [assly_.no] . 

dol [lot_no] . units_lateness_of_lot: = 
(current_supply_period - current_demand_period) * 
current_demand_qty; 
total_units_lateness: = 
total_units_lateness + 

(current_supply_period - current_demand_period) * 
current_demand_qty; 
end; (* then * ) 

current_supply_qty: =current_supply_qty-current_demand_qty; 
end; (* while current_supply_qty > 0 *) 
end; (* for lot_no *) 
end; (* for assly_no *) 

(* end of calculating units_lateness *) 

(* calculating me utilisation statistics *) 

for machine_.no : =1 to max_mc_no do 
begin (* for machine_no *) 
agm [machine_.no] . cu: = 
agm [machine_.no] . tpt/ 

(clock_time - collect_statistics_from_time) ; 
end; (* for machine_no *) 

for machine_.no : =16 to 20 do 
begin (* for machine_no *) 
aam [machine_no] . cu: = 
aam [machine_.no ] . tpt/ 

(clock_time ~ collect_statistics_from_time) ; 
end; (* for machine_no *) 

(* end of calculating me utilisation statistics *) 
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if (total_units_lateness < best_sol) then 
begin 

rewrite (potty_output, ’ s. res’ ) ; 

writeln(potty_output total aggregate lateness =’ , 

total_aggregate_lateness) ; 

writelnlpot ty output , total se tup cost — ’ , total setup cost ) * 

writeln(potty_output , ’ total units lateness =’ , total_units_lateness) ; 
writeln(potty_output, ’ total inventory carrying cost =’ , 

total_inventory_carrying_cost) ; 

writeln(potty_output, ’ theoretical inventory carrying cost =’ , 

theoretical_inventory_carrying_cost ) ; 

sum_ut i 1 : =0 . 0 ; 

for general_mc: = 1 to max_no_of_proces_on_any_assly do 

sum_util: =sum_util+agm[general_mc] . cu; 

av_util : =sum_util/max_no_of_proces_on_any_assly; 

for general_mc:= 1 to max_no_of_proces_on_any_assly do 

writeln(potty_output , 

’Utilisation for general machine’ ,general_mc: 2, ’ =’ , 

agm[general_mc] . cu: 5) ; 

writeln(potty_output, ’ Average Utilisation is’,’ =’ , av_util : 5) ; 

surnjrange : =0 ; 

for mc_no : = 16 to 20 do 

begin 

range [mc_no] : =latest_lot_time [mc_no] - 

earliest_lot_time [mc_no] ; 

writeln(potty_output, ’ coordination_range for assly centre’ ,mc_no: 2, 
’is’, range [ mc_no ] ) ; 
sum_range: =sum_range+range [mc_no] ; 

end; 

coordination_range: =sum_range/ (20-15) ; 

if first_run then 

prev_range: =coordination_range; 

writeln(potty_output, ’ previous coordination range is’,’ =’ , 

prev_range : 5 ) ; 

writeln(potty_out.put ,’ improved coordination range is’,’ =’ , 

coordination_range :5); 

close (potty_output ) ; 
rewrite(s_dat , ’ s. dat’ ) ; 

writeln(s_dat , total_aggregate_lateness) ; 
writeln(s_dat , total_setup_cost ) ; 
writeln(s_dat , total_units_lateness) ; 
writeln(s_dat , total_inventory_carrying_cost) ; 
writeln(s_dat , theoretical_inventory_carrying_cost) ; 

writeln(s_dat , av_util ) ; 
writeln(s_dat , coordination_range ) ; 
close (s_dat) ; 

for general_mc:= 1 to max_no_of_proces_on_any_assly do 
writeln( ’ Utilisation for general machine’ ,general_mc: 2, ’ =’ , 

agm[general_mc] . cu: 5) ; 

writelnC ’ Average Utilisation is’,’ =’ ,av_util:5) ; 

end; 

end; (* collect_and_print_statistics *) 
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procedure set_parameters; 
var 

assly, 
assly_no, 
lot_no, 
no_of_lots , 

mc_no : integer; 

begin 

for assly_no:= 1 to max_assly_no do 
if abda[assly_no] . immediately_needs_processing then 

begin 

assly: =abda [ ass ly_no ] . assly_no_af ter processing; 

critical_child [assly] : =assly_no; 

no_of_lots : =array_f or_output_of_w_w [assly_no ] . 

no_of_orders_placed; 
for lot_no:= 1 to no_of_lots do 
begin 

time_of_reaching [assly_.no , lot_no] :=0; 

end; 

coordinat.ion_index[assly_no] : =1 ; 
for mc_no:= 16 to 20 do 
latest_assly_time [assly_no, mc_no] : =0; 
end; 

for mc_no:= 16 to 20 do 
begin 

latest_lot_time [mc_no] : =0; 
earliest_lot_time [mc_no] : =0; 
max_t ime [ mc_no ] : =0 ; 
assly_count. [mc_no] : =0; 
presen t_mean_t ime [mc_no] : =0; 

end; 

end; 

procedure next_set_parameters; 
var 

assly, 
assly_.no, 
lot_no, 
no_of_lots , 

mc_no : integer; 

begin 

for assly_no:= 1 to max_assly_no do 

if abda [assly_no] . immediately_needs_processing then 

begin 

no_of_lots: =array_for_output_of_w_w[assly_no] . 

no_of _orders_placed ; 

for lot_no:= 1 to no_of_lots do 
begin 

t ime_of peaching [ ass ly_no , lot_.no] :=0; 

end; 

for mc_no:= 16 to 20 do 
latest_assly_time [assly_no, mc_no] : =0; 
end; 



for mc_no:= 16 to 20 do 
begin 

latest_lot_time [mc_no] :=0; 
earliest_lot_time [mc_no] : =0; 
max_t ime [mcjno ] : =0 ; 

prev_mean_time [mc_no] : =present_mean_time [me no] : 
pr esent_mean_t ime [ mcjno ] : =0 ; 

end; 

end; 


procedure run; 
var 

a, b: integer; 
begin 

read_input_data ; 
do_explosion; 
do_lot_sizing; 
initialisations; 
best_sol : =inf ini ty ; 
set_parameters ; 
prepare_queue_of_lots ; 

prepare_event.s_regarding_entry_of_components; 
prepare_collect_statistic_event ; 

(* decide_times_of_entry_in_t.he_shop; *) 
open_print_the_event_list ; 

while (clock_t.ime < simulation_time) and (flag = false) do 
begin 

f ind_the_next._executable_event ; 

if (flag=false) then 

begin 

pr int_the_event_l ist ; 
execute_the_event ; 
end; 
end; 

close_print_the_event_l ist ; 
collect_and_print_statistics; 

end; 


procedure next_run; 
var 

a,b: integer; 
begin 

do_lot_sizing; 
initialisations ; 
next_setjparameters; 
prepare_queue_of_lots ; 

prepare_events_regarding_entry_of_components; 

prepare_collect_statistic_event; 

open_pr int_the_event_l ist ; 

while (clock_time < simulation_time) and (flag = false) do 
begin 

f ind_the_next_executable_event ; 

if (flag=false) then 

begin 
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print_the_event_list ; 
execute_the__event ; 
end; 
end; 

closejprint_the__event_list ; 
col lec t_andjprint_s tat i sties; 

end; 


begin (* main *) 
times: =0; 
f irst__run: =true; 
run; 

f irst_run: =false; 

(* while ( total_units_lateness <= best_sol) and (times < 10) do*) 

while ( times<20) do 

begin 

times: =t imes+1 ; 

writeln( times, * iteration’ ) ; 

if (total_units_lateness<best_sol) then 

begin 

writelnt’ improved’ ) ; 

best._sol : =total_units_ lateness ; 

end; 

next_run; 

end; 

end. 




TABLE B.l ROUTING DETAIL OF COMPONENTS 
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Assly No. 

No. of demands 

Period 

1 

16 

1 ,4,7,1 2,1 5,1 7,20,22,26,30,35,39,42,46,50 


TABLE B.2 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

705149 

748046 

692481 

671361 

TOtAL 

SETUP 

COST 

1645 

1721 

1645 

1645 

TOTAL 

UNITS 

LATENESS 

756 

712 

756 

701 

TOTAL 

INVENTORY 

CARRYING 

COST 

24400400 

27400500 

24400400 

22600400 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

1 1 1 07.7 

10528 

11107.7 

11107.7 

AVERAGE 

UTILISATION 

0.7443 

0.757817 

0.7443 

0.774761 

AVERAGE 

COORDINATION 

RANGE 

713.5 

371 .25 

713.5 

283.6 


TABLE B.3 


RESULTS (PROBLEM 1) 




















Assly No. 

No. of demands 

Period 

1 

9 

1 ,4,7,1 9,29.35,41 ,47.50. 

3 


7,14,21,35,42 

4 

5 

6,18,30,36,48 

5 

5 

9,18,27,36,45 


TABLE B.5 DEMAND DETAILS 


” PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

778019 

791597 

782741 

750593 

TOTAL 

SETUP 

COST 

740 

780 

740 

740 

TOTAL 

UNITS 

LATENESS 

426 

345 

394 

312 

TOTAL 

INVENTORY 

CARRYING 

COST 

23800000 

26400000 

23800000 

23400000 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

“1434^ 

11 434.3 

11434.3 

11434.3 

AVERAGE 

UTILISATION 

0.673021 

0.780576 

0.692023 

0.731424 

AVERAGE 

COORDINATION 

RANGE 

479 

312 

279 

184 


TABLE B.6 


RESULTS (PROBLEM 2) 
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Assly No. 

No. of demands 

Period 

1 

12 

1 ,3,7,13,17.27.31 .33,37,41 .48.50 

3 

4 

6,16,36,46 

4 

5 

7,13,20,27,37 

9 

5 

2,11,20,29,39 


TABLE B.8 DEMAND DETAILS 


TOTAL 

AGGREGATE 

LATENESS 

995432 

995540 

1020560 

1021470 

TOTAL 

SETUP 

COST 

1265 

1265 

1265 

1265 

TOTAL 

UNITS 

LATENESS 

802 

802 

744 

698 

TOTAL 

INVENTORY 

CARRYING 

COST 

44199900 

44199900 

44199900 

41999900 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

15393* 

15393 

15393 

15393 

AVERAGE 

UTILISATION 

0.713262 

0.713262 

0.738163 

0.816929 

AVERAGE 

COORDINATION 

RANGE 

746.25 

453.5 

746.25 

307.6 


TABLE B.9 


RESULTS (PROBLEM 3) 
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Assly No. 

No. of demands 

Period 

1 

16 

1.4.7, 1 0,12,1 5,17,20,22,26.30.35,39,42,46.50 


TABLE B.1 1 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

812520 

746240 

811881 

712532 

TOTAL 

SETUP 

COST 

1461.4 

1628 

1461.4 

1461.4 

TOTAL 

UNITS 

UTENESS 

1288 

1650 

1288 

1134 

TOTAL 

INVENTORY 

CARRYING 

COST 

37200200 

37000600 

37200200 

33400200 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

8737.58 

7249.93 

8737.58 

8737.58 

AVERAGE 

UTILISATION 

AVERAGE 

COORDINATION 

RANGE 

0.753628 

805.975 

0.767539 

0.742674 

0.753628 

570,175 

810.975 

553.28 


TABLE B.12 


RESULTS (PROBLEM 4) 




FIGURE 13.5 PROBLEM 5 


Component No 

Routing Sequence 

2 

1 4 3 4 5 

5 

1 4 2 4 3 

6 

3 4 5 

11 

1 4 2 

12 

3 4 4 

13 

2 4 1 

14 

4 4 3 

17 

3 -> 4 

18 

1 4 2 

21 

5 4 3 4 2 

22 

4 4 1 


TABLE B. 13 ROUTING DETAIL OF 
COMPONENTS 
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Assly No. 

1 

No. of demands 

13 

— Period 1 

_ _. _ 1 > 4 . 6 . 12, 15, 16, 19, 26.31 37 4R4R«;n 

12 

4 

2,18,26,43 

17 

5 

7,14,21 ,35,42 


TABLE B.14 DEMAND DETAILS 


PARAMETER 

TOTAL 

AGGREGATE 

LATENESS 

ORIGINAL 

hT 

H2 

H3 

768486 

704840 

750229 

668486 

TOTAL 

SETUP 

COST 

1548 

1995 

1548 

1548 

TOTAL 

UNITS 

LATENESS 

1489 

2375 

1558 

1489 

TOTAL 

INVENTORY 

CARRYING 

COST 

35000300 

52800700 

33800300 

33000300 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

8737.58 

7249.93 

8737.58 

8737.58 

AVERAGE 

UTILISATION 

0.786595 

0.817647 

0.768112 

0.786595 

“AVERAGE 

COORDINATION 

RANGE 

910'5 

947.25 

934.25 

620.8 


TABLE B.A5 RESULTS (PROBLEM 5) 



1 



TABLE B. 16 ROUTING DETAIL OF COMPONENTS 


Demand details 
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Assly No. 

1 

No. of demands 

13 

— Period — i 

— 1.4,7,10,12,1 6,19,26.31.35,41 ,47 50 

4 

4 

1.21,31,41 

6 

5 

7,22.32,37,42 


TABLE B.17 DEMAND DETAILS 


““PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

954103 

958472 

978046 

972650 

TOTAL 

SETUP 

COST 

2132 

2160 

2132 

2132 

TOTAL 

UNITS 

LATENESS 

1917 

1762 

1762 

1337 

TOTAL 

INVENTORY 

CARRYING 

COST 

56200400 

53800500 

58200400 

76400400 

THEORETICAL"' 

INVENTORY 

CARRYING 

COST 

8804.97 

8804.97 

8804.97 

88049700 

" AVERAGE 
UTILISATION 

0.849468 

0.887055 

0.870687 

0.811669 

AVERAGE 

COORDINATION 

RANGE 

888 

492.5 

890.25 

488.4 


TABLE B.18 


RESULTS (PROBLEM 6) 
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Assly No. 

No. of demands 

Period 1 

1 

14 ^ 

j 

- 1 .7,9,14,1 6,21 ,23.27.31,39.41.47.50 

3 

4 

1,21,41,49 

5 

4 

2,22,42,47 

15 

4 

3,11,25,31 

23 

4 

7,34,36,47 


TABLE B.20 DEMAND DETAILS 


PARAMETER 

ORIGINAL, 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

1365700 

1327340 

1365700 

1265000 

TOTAL 

SETUP 

COST 

1932 

4 

2006 

1932 

1932 

TOTAL 

UNITS 

LATENESS 

"“1392 

796 

1392 

1317 

TOTAL 

INVENTORY 

CARRYING 

COST 

~ 54800300 

74600100 

54800300 

66000100 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

8670.92 

8536.05 

8670.92 

8670.92 

AVERAGE 

UTILISATION 

0.64G61 7 

0.639394' 

0.646617 

0.655988 

AVERAGE 

COORDINATION 

RANGE 

810.75 

565.25“ 

810.75 

466.4 


TABLE B.21 


RESULTS (PROBLEM 7) 
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Assly No. 

No. of demands 

Period 

1 

12 

1 ,4,7,1 1 ,1 5,1 9,24,28,32,38,46,50 

4 

5 

3,7,17,23,43 

5 

5 

2,7,25,39,46 


TABLE B.23 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

1049110 

1049270 

1054680 

992817 

TOTAL 

SETUP 

COST 

1804 

1804 

1804 

1804 

TOTAL 

UNITS 

LATENESS 

772 

772 

728 

684 

TOTAL 

INVENTORY 

CARRYING 

COST 

52600400 

52600400 

53600400 

59000300 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

8873.66 

8873.66 

8873.66 

8873.66 

AVERAGE 

UTILISATION 

0.676471 

0.676471 

0.680639 

0.645474 

AVERAGE 

COORDINATION 

RANGE 

556 

548 

557.5 

161.6 


TABLE B.24 


RESULTS (PROBLEM 8) 






FIGURE B. 9 PROBLEM 9 


Components 

Routing Sequence 

2 

1 -» 4 -> 7 

5 

1 2 3 4 

6 

5 6 7 -> 4 

12 

1 -> 2 -» 3 

13 

2 -> 3 -» 4 

14 

7 -4 4 -» 2 

15 

3 4 4 -> 5 

16 

4 -> 5 -> 6 

19 

5 -» 6 -> 7 

20 

6 -> 7 -> 1 


TABLE B. 25 ROUTING DETAIL OF COMPONENTS 
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Assl>^ n d 

No. of demands 

Period 


10 

1,4,11,19,22,28,33,39,44,50 

59 

6 

4,8,14,18,21,29 


TABLE B _ 26 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

total 

aggregate 

LATENESS 

936183 

986948 

975974 

913242 

total 

SETUP 

COST 

1548 

1683 

1548 

1548 

TOTAL 

UNITS 

LATENESS 

1260 

1118 

976 

976 

TOTAL 

INVENTORY 

CARRYING 

COST 

34200100 

45800100 

37400100 

31100000 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

15393 

15393 

15393 

15393 

AVERAGE 

UTILISATION 

0.55773 

0.586314 

0.587972 

0.560045 

AVERAGE 

COORDINATION 

RANGE 

835 

571 .25 

835 

274.8 


TABLE B.27 


RESULTS (PROBLEM 9) 





TABLE B. 28 ROUTING DETAIL OF COMPONENTS 
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Assly No. _ 

No, of demands 

Period 



1 

13 

1 .3.7,1 1 ,16,21 ,24,27,31 ,35,41 ,46,50 

7 

17.27,37,41,43,47.50 

9 .... J 



4 n 

4 

3.11,29,31 

YC. 

L Ji- - 

4 

1,21,29,43 


TABLE B.29 DEMAND DETAILS 



HI 

748046 

H2 

H3 

ORIGINAL 

705149 

692481 

671361 



1*791 

1645 

1645 

1645 

I / C 1 




“Tl'2" 

756 

^toT 

756 




r\~T A AACnn 

“24400400 

22600400 

24400400 

274UUOUU 

- “THorT 

11107.7 

10528 

11107.7 


. 0.7443 

' " 0.757817 

0.7443 

0.774761 

- — — — “ “ ?R3ir 

713.5 

37T.2S 

L - 

_J 


TABLE B.30 RESULTS (PROBLEM 10) 



FIGURE B. 11 PROBLEM 11 



TABLE B. 31 ROUTING DETAIL OF COMPONENTS 












Assly No. 

No.of demands 

Period 

1 

9 

1,4,9,17,21,28,39.44.50 

3 

4 

9,19,29,39 

4 

6 

3,15,22,36,42,46 

17 

6 

8,18,24,30,39,46 


TABLE B.32 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

1210380 

1210380 

1210380 

1191560 

TOTAL 

SETUP 

COST 

1538 

1538 

1538 

1538 

TOTAL 

UNITS 

LATENESS 

382 

382 

382 

349 

TOTAL 

INVENTORY 

CARRYING 

COST 

48800300 

48800300 

48800300 

50600300 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

... 

8354.61 

8354.61 

8354.61 

AVERAGE 

UTILISATION 

0.765354 

0.765354 

0.765354 

0.708108 

AVERAGE 

COORDINATION 

RANGE 

482.5 

74.75 

482.5 

69.6 


TABLE B.33 


RESULTS (PROBLEM 11) 



























Assly No. 

No.of demands 

Period 

1 

10 

1,3,7.11,17,19,25.35.45.50 

6 

5 

9,19,27,35,50 

8 

5 

2,4,32,37,45 

13 

^ 4 

4,21,37,46 


TABLE B.35 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

871848 

871848 

871848 

774585 

TOTAL 

SETUP 

COST 

1138 

1138 

1138 

1138 

TOTAL 

UNITS 

LATENESS 

" !1A •••" ' 11 " " ' ' IU ' 

1474 

1474 

1474 

1282 

TOTAL 

INVENTORY 

CARRYING 

COST 

23200400 

23200400 

23200400 

23600300 

“THEORETICAL 

INVENTORY 

CARRYING 

COST 

16233.8 

16233.8 

16233.8 

16233.8 

AVERAGE ~ 

UTILISATION 

0.711028 

0.711028 

0.711028 

0.783936 

AVERAGE 

COORDINATION 

RANGE 

694" 

142.25 

694 

87.2 


TABLE B.36 


RESULTS (PROBLEM 12) 




Component No. 

Routing Sequence 

2 

1 -)4-)5 

6 

2 -» 3 

7 

1 4 

8 

— 5 4 

14 

5 4 1 

15 

3 5 

16 

. 4 -> 2 

17 

— 3-74 

18 

“ 1 -»"3 

21 

~ m 

22 

2^4 


TABLE B. 37 ROUTING DETAIL OF COMPONENTS 


Demand details 
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Assly No. 

No.of demands 

Period 

1 

10 

6,11,17,22.27,31.38,41,46,50 

4 

6 

7,11,19,27,31,41 

8 

6 

4,14,19,27,35,46 

10 

6 

3,13,27,37,41,50 


TABLE B.38 DEMAND DEATAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

921864 

904073 

912091 

782381 

TOTAL 

SETUP 

COST 

1392 

1412 

1392 

1392 

TOTAL 

UNITS 

LATENESS 

1937 

1462 

1651 

1418 

TOTAL 

INVENTORY 

CARRYING 

COST 

211000 

242900 

92763 

8.5643 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

8372.63 

8372.63 

8372.63 

8372.63 

AVERAGE 

UTILISATION 

0.814752 

0.793945 

0.803269 

0.827234 

AVERAGE 

COORDINATION 

RANGE 

1142 

279.5 

1 1 67.5 

136.8 


TABLE B.39 


RESULTS (PROBLEM 13) 





FIGURE B. 14 PROBLEM 14 


Component No. 

Routing Sequence 

2 

5 -» 3 -» 1 

7 

1 -> 2 -> 3 

8 

2 4 

9 

2 -> 3 4 

10 

3 4 

15 

3 -> 4 

16 

4 -» 5 

17 

1 2 

18 

3 4 

21 

4 -» 5 -> 1 

22 

5 -» 3 -> 4 


TABLE 40 ROUTING DETAIL OF COMPONENTS 
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Assly No. 

No. of demands 

Period 

1 

10 

1,9,21,26,29,31,38,40,46,50 

4 

7 

2,9,21,29,33,41,47 

6 

5 

6,16,26,36,46 

14 

5 

9,14,22,29,38 


TABLE B.41 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

938586 

938586 

943955 

804328 

TOTAL 

SETUP 

COST 

1472 

1472 

1472 

1472 

TOTAL 

UNITS 

LATENESS 

1524 

1524 

1394 

1329 

TOTAL 

INVENTORY 

CARRYING 

COST 

24600400 

24600400 

24400400 

25400400 

! 

i 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

9548.41 

9548.41 

9548.41 

9548.41 

AVERAGE 

UTILISATION 

0.681069 

0.681069 

0.741236 

0.743636 

AVERAGE 

COORDINATION 

RANGE 

1027.5 

661.5 

1030.25 

513.2 


TABLE B.42 


RESULTS (PROBLEM 14) 







FIGURE B. 15 PROBLEM 15 


Component No. 

Routing Sequence 

2 

2 -> 3 -> 5 

6 

5 -> 3 2 

7 

3 -> 4 5 

8 

1 2 3 

13 

1 2 3 4 

14 

2 3 4 5 

15 

5 -» 4 -> 3 -> 2 

16 

3 4 -> 1 5 

20 

1 -> 2 

21 

3 5 

22 

4^2 


TABLE B. 43 ROUTING DETAIL OF COMPONENTS 























Assly No. 

No. of demands 

Period 

1 

9 

3,9,15,21,29,32,39,41,48 

3 

3 

9,29,37 

5 

6 

2,11,17,24,31,39 

12 

5 

5,15,25,35,45 


TABLE B.44 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

880565 

735010 

815694 

642816 

TOTAL 

SETUP 

COST 

1894 

1870 

1894 

1894 

TOTAL 

UNITS 

LATENESS 

921 

2171 

1311 

270 

TOTAL 

INVENTORY 

CARRYING 

COST 

311000 

337000 

237000 

213220 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

8170.46 

7146.09 

8170.46 

: 

8170.46 

AVERAGE 

UTILISATION 

0.815569 

0.822999 

0.726952 

0.851916 

AVERAGE 

COORDINATION 

RANGE 

858.5 

907 

928.5 

300.2 


TABLE B.45 


RESULTS (PROBLEM 15) 







FIGURE B. 16 PROBLEM 16 


Component No. 

Routing Sequence 

2 

2 -> 4 6 

7 

5 -> 2 -> 3 

8 

1 ^ 4 -> 2 

9 

3 -> 4 -» 5 

10 

3 -> 4 ^ i 

15 

1 2 

16 

3 4 

17 

5 6 

18 

•2 -> 4 


TABLE B. 46 ROUTING DETAIL OF COMPONENTS 
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Assly No. 

No. of demands 

Period 

1 

10 

1,3,7,17,21,26,29,31,39,47 

5 

4 

9,27,33,47 


TABLE B.47 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

583868 

644377 

587148 

560322 

TOTAL 

SETUP 

COST 

895 

1288 

895 

895 

TOTAL 

UNITS 

LATENESS 

582 

690 

654 

548 

TOTAL 

INVENTORY 

CARRYING 

COST 

21200100 

26600200 

21200100 

21200000 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

8493.47 

10186.9 

8493.47 

8493.47 

AVERAGE 

UTILISATION 

0.748428 

0.745731 

0.638027 

0.802671 

AVERAGE 

COORDINATION 

RANGE 

299.75 

166.5 

318.25 

86.6 


TABLE B.48 


RESULTS (PROBLEM 16) 






FIGURE B. 17 PROBLEM 17 


Component No. 

Routing Sequence 

2 

1 -» 3 -» 6 

7 

1 3 -> 5 

8 

2 4 -> 6 

9 

6 4 2 

10 

5 3 1 

15 

1 -» 2 

16 

3 -> 4 

17 

5^6 

18 

2 -> 4 -> 6 


TABLE B. 49 ROUTING DETAIL OF COMPONENTS 



















Assly No. 

No. of demands 

Period 

1 

11 

1 .3,6,1 1 ,1 5,21 ,27,31 ,35,42,50 

4 

4 

4,14,34,44 

6 

5 

6,16,26,36,46 


TABLE B.50 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

537844 

584176 

543596 

521826 

TOTAL 

SETUP 

COST 

1372 

1781 

1372 

1372 

TOTAL 

UNITS 

LATENESS 

810 

696 

753 

696 

TOTAL 

INVENTORY 

CARRYING 

COST 

22200100 

23000100 

22200100 

21400100 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

4254.23 

4254.23 

4254.23 

4254.23 

AVERAGE 

UTILISATION 

0.623706 

0.613054 

0.680936 

0.738272 

AVERAGE 

COORDINATION 

RANGE 

368.75 

247.75 

388.25 

157.2 


TABLE B.51 


RESULTS (PROBLEM 17) 





FIGURE B. 18 PROBLEM 18 


Component No. 

Routing Sequence 

2 

2 -> 4 5 -> 6 

6 

6 5 4 

7 

3 -> 5 -> 4 

8 

4 -» 5 3 

14 

2 1 

15 

5 -> 2 

16 

4-^6 

17 

3 -> 6 

18 

1 5 


TABLE B. 52 ROUTING DETAIL OF COMPONENTS 
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Assly No. 

No. of demands 

Period 

1 

12 

1 ,3,7,1 1 ,1 9,21 ,28,32,36,42,44,50 

4 

6 

9,14,21,27,35,45 

5 

5 

3,9,16,36,45 


TABLE B.53 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

742095 

747484 

741962 

727928 

TOTAL , 

SETUP 

COST 

1275 

1507 

1275 

1275 

TOTAL 

UNITS 

LATENESS 

600 

634 

600 

532 

TOTAL 

INVENTORY 

CARRYING 

COST 

21000400 

37000300 

21000400 

35200200 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

4816.42 

3934.46 

4816.42 

4816.42 

AVERAGE 

UTILISATION 

0.701754 

0.623341 

0.701754 

0.685289 

AVERAGE 

COORDINATION 

RANGE 

432.5 

220.25 

432.5 

142.8 


TABLE B.54 


RESULTS (PROBLEM 18) 





FIGURE B. 19 PROBLEM 19 


Component No. 

Routing Sequence 

2 

1 -> 3 -» 5 -> 7 

6 

1 -> 2 3 

7 

4 -> 5 6 

8 

7 -> 2 -» 1 

14 

7 -> 3 -> 4 

15 

4 5 6 

16 

3 -> 4 -> 1 

17 

7 -> 6 -» 5 

18 

6 4 -» 2 

23 

3 1 -» 2 

24 

5 -» 6 -» 7 

25 

4 -> 6 -> 1 

26 

2 -> 5 4 3 


TABLE B. 55 ROUTING DETAIL OF COMPONENTS 




























Assly No. 

No. of demands 

Period 

1 

10 

1,4,9,17,21,28,33,40,46,50 

4 

9 

3,9,13,19,23,29,33,39,46 

5 

5 

5,15,25,35,45 

ii 

3 

21 ,31,41 


TABLE B.56 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

HI 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

1690930 

1758050 

1705070 

1 640580 

TOTAL 

SETUP 

COST 

2534 

3918 

2534 

2534 

TOTAL 

UNITS 

LATENESS 

1283 

1632 

1283 

973 

TOTAL 

INVENTORY 

CARRYING 

COST 

53600400 

65600600 

51400400 

53200100 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

13838.3 

13838.3 

13838.3 

13838.3 

AVERAGE 

UTILISATION 

0.759847 

0.771538 

0.743453 

0.840568 

AVERAGE 

COORDINATION 

RANGE 

876.75 

559.25 

889.5 

367.2 


TABLE B.57 


RESULTS (PROBLEM 19) 






FIGURE B. 20 PROBLEM 20 


Component No. 

Routing Sequence 

2 

2 4 6 

5 

7 -> 6 -> 5 

6 

4 3 2 

15 

1 2 

16 

3 -> 4 

17 

5 -> 6 

18 

7 1 

19 

3 2 

20 

4 5 

21 

6 -> 1 

22 

2 -s. 4 

25 

6 -> 3 

26 

7 4 


TABLE B. 58 ROUTING DETAIL OF COMPONENTS 

























Assly No 

No. of demands 

Period 

1 

12 

1 ,5,9,1 3,1 7,21 ,27,31 ,35,39,43,47 

3 

9 

6,9,11,17,23,29,33,39,47 

4 

5 

3,13,23,33,43 

12 

3 

19,38,47 


TABLE B.59 DEMAND DETAILS 


PARAMETER 

ORIGINAL 

Hi 

H2 

H3 

TOTAL 

AGGREGATE 

LATENESS 

1044140 

1120730 

1046690 

1023130 

TOTAL 

SETUP 

COST 

1833 

2217 

1833 

1833 

TOTAL 

UNITS 

LATENESS 

1214 

1159 

1214 

994 

TOTAL 

INVENTORY 

CARRYING 

COST 

38200900 

44201300 

38200900 

34600600 

THEORETICAL 

INVENTORY 

CARRYING 

COST 

22681 .6 

22117 

22681 .6 

22681 .6 

AVERAGE 

UTILISATION 

0.728353 

0.818681 

0.727625 

0.855085 

AVERAGE 

COORDINATION 

RANGE 

747.75 

523.25 

760 

317.6 


TABLE B.60 


RESULTS (PROBLEM 20) 




IMPROVE 

PERCENTAGE 

MENT OVER ORIGINAL 

HI 

H2 

H3 

6.401055 

33.2876 

34.70185 

19.01408 

7.511737 

26.76056 

0 

7.23192 

12.96758 

-28.1056 

0 

11.95652 

-59.503 

-4.63398 

0 

8.08555 

8.08555 

30.25561 

42.81609 

0 

5.387931 

0 

5.699482 

11.39896 

11.26984 

22.53968 

22.53968 

5.820106 

0 

7.275132 

0 

0 

8.638743 

0 

0 

13.02578 

24.52246 

14.7651 

26.79401 

0 

8.530184 

12.79528 

-135.722 

-42.3453 

70.68404 

-18.5567 

-12.3711 

5.841924 

14.07407 

7.037037 

14.07407 

-5.66667 

0 

11.33333 

-27.2019 

0 

24.16212 

4.530478 

0 

18.12191 

-7.21111 

2.766895 

18.43575 


TABLE B.61 PERCENTAGE IMPROVEMENT (TOTAL UNITS LATENESS) 


L 

HI Vs H2 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

m 

value from table 

0.95 

38 

1.15128 

1.686 


| H2 Vs H3 | 

k 

degrees of 

FREEDOM 

|T| 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.995 

38 

3.41075 

2.713 



HI Vs H3 

k 

DEGREES OF 
FREEDOM 

|T| 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.995 

38 

2.93549 

2.713 


TABLE B.62 t TEST DATA (TOTAL UNITS LATENESS) 








PERCENTAGE IMPROVEMENT 1 

I OVER ORIGINA 


HI 

H2 

H3 

-0.04794 

23.87482 

26.41635 

-1.7452 

-0.60693 

3.525107 

-0.01085 

-2.52433 

-2.61575 

8.157338 

0.078644 

12.30591 

8.281999 

2.37571 

13.0126 

-0.45792 

-2.50948 

-1.94392 

2.808816 

0 

7.373508 

-0.01525 

-0.53093 

5.365786 

-5.42255 

-4.25034 

2.450482 

-6.0834 

1.7965 

4.791611 

0 

0 

1.554884 

0 

0 

11.15596 

1.929894 

1.060135 

15.13054 

0 

-0.57203 

14.30428 

16.52973 

7.366975 

26.9996 

-10.3635 

-0.56177 

4.032761 

-8.61439 

-1.06946 

2.978187 

-0.72619 

0.017922 

1 .909055 

-3.96941 

-0.83623 

2.977651 

-7.33522 

-0.24422 

2.012182 

-0.3542 

1.14325 

7.686839 


TABLE B.63 PERCENTAGE IMPROVEMENT (TOTAL AGGREGATE LATENESS) 


HI Vs H2 | 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.95 

38 

0.80593 

1.686 


H2 Vs H3 | 

k 

DEGREES OF 
FREEDOM 

|T| 

CALCULATED VALUE 

m 

VALUE FROM TABLE 

0.995 

38 

2.96734 

2.713 


HI Vs H3 I 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.995 

38 

4.32772 

2.713 


TABLE B.64 t TEST DATA (TOTAL AGGREGATE LATENESS) 







PERCENTAGE 
IMPROVEMENT OVER 
ORIGINAL 

HI 

H2 

H3 

- 11.4507 

-29.0037 

31.35998 

— 10.9244 

0 

1 .680672 

0 

0 

4.977387 

0.536556 

0 

10.215 

-50.8578 

3.428542 

2.366843 

4.270254 

-3.55869 

-31.2713 

-36.1308 

0 

-20.4375 

0 

-1.90113 

-10.0744 

-33.918 

-9.3567 

16.84514 

-12.2953 

0 

7.376928 

0 

0 

-3.6885 

0 

0 

-1.72368 

-15.1185 

56.03649 

99.99077 

0 

0.812995 

-4.09829 

-8.36013 

23.79421 

10.03376 

-25.472 

0 

0.000472 

-3.60359 

0 

3.603587 

-76.1885 

0 

-67.6168 

-22.3883 

4.104447 

-3.50133 

-15.7075 

0 

9.424647 

-15.8804 

2.217821 

2.773173 


TABLE B.65 PERCENTAGE IMPROVEMENT (TOTAL INVENTORY CARRYING COST) 


I HI Vs H2 1 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.995 

38 

3.23572 

2.713 


. . _ 

H2 Vs H3 1 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

t(k> 

VALUE FROM TABLE 

0.95 ! 

38 

0.07454 

1.686 


H?VsH3 I 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.995 

38 

3.33501 

2.713 


TABLE B.66 t TEST DATA (TOTAL INVENTORY CARRYING COST) 






PERCENTAGE 
IMPROVEMENT OVER 
ORIGINAL 

HI 

H2 

H3 

0.587775 

39.44973 

32.09222 

15.98093 

2.823389 

8.677738 

0 

3.491144 

14.53421 

1 .845871 

-1.4535 

0 

3.947648 

-2.34975 

0 

4.424769 

2.497916 

-4.44973 

-1.11704 

0 

1.449235 

0 

0.616139 

-4.58216 

5.125061 

5.422337 

0.415075 

1.816069 

0 

4,09257 

0 

0 

-7.47968 

0 

0 

10.25389 

-2.55378 

-1.40939 

1.532 

0 

8.8342 

9.186588 

0.91102 

-10.8657 

4.456643 

-0.36036 

-14.7511 

7.247591 

-1.70786 

9.175798 

18.36859 

-11.1739 

0 

-2.34626 

1.538599 

-2.15754 

10.62332 

12.40168 

-0.09995 

17.3998 

1.583326 

1.961191 

6.073582 


TABLE B.67 PERCENTAGE IMPROVEMENT (AVERAGE UTILIZATION) 


[ — 

HI Vs H2 | 

k 

DEGREES OF 
FREEDOM 

|T| 

CALCULATED VALUE 

Kk) 

VALUE FROM TABLE 

0.95 

38 

0.46742 

1.686 


H2 Vs H3 

k 

DEGREES OF 
FREEDOM 

|T| 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.9995 

38 

4.23881 

3.57 



HI Vs H3 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.9995 

38 

5.55442 

3.57 


TABLE B.68 t TEST DATA (AVERAGE UTILIZATION) 







PERCENTAGE IMPROVEMENT 
OVER ORIGINAL 

HI 

H2 

H3 

0.349502 

-0.86316 

73.65812 

34.8643 

41 .75365 

61.58664 

39.22948 

0 

58.78057 

29.25649 

-0.62037 

31 .35271 

-4.03624 

-2.60846 

31.81768 

44.53829 

-0.25338 

45 

30.2806 

0 

42.47302 

1.438849 

-0.26978 

70.93525 

31.58683 

0 

67.08982 

47.96776 

0 

60.25228 

84.50777 

0 

85.57513 

79.50288 

0 

87.43516 

75.52539 

-2.23292 

88.02102 

35.62044 

-0.26764 

50.05353 

-5.64939 

-8.15376 

65.03203 

44.45371 

-6.17181 

71.10926 

32.81356 

-5.28814 

57.36949 

49.07514 

0 

66.98266 

36.21329 

-1.45423 

58.11805 

30.0234 

-1.63825 

57.52591 

35.8781 

0.596588 

61.50842 


TABLE B.69 PERCENTAGE IMPROVEMENT 

(AVERAGE COORDINATION RANGE) 


HI Vs H2 | 

k 

DEGREES OF 
FREEDOM 

m 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.9995 

38 

5.789717 

3.57 


_ _ _____ 

H2 Vs H3 | 

k 

DEGREES OF 
FREEDOM 

|T| 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.9995 

38 

14.3427 

3.57 


| HI Vs H3 

k 

DEGREES OF 
FREEDOM 

|T| 

CALCULATED VALUE 

t(k) 

VALUE FROM TABLE 

0.9995 

38 

4.20595 

3.57 


TABLE B.70 t TEST DATA (AVRAGE COORDINATION RANGE) 







